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1.0  INTRODUCTION 

The  effects  of  materia!  outgassing  are  of  continuing  concern  to  satellite  designers  since 
they  are  being  asked  to  extend  the  usable  lifetime  of  space  satellites.  The  lifetime  is  most 
often  directly  related  to  contamination  of  spacecraft  critical  surfaces.  During  the  past  few 
years  considerable  progress  has  been  made  at  AEDC  in  the  area  of  contamination  and 
optical  property  measurements.  Both  experimental  and  analytical  efforts  have  been  made  in 
the  determination  of  optical  properties  used  in  characterizing  contaminants. 

The  basic  property  required  for  determining  optical  contamination  effects  is  the  complex 
refractive  index  n(?)  =  n(p)  +  ik  (p),  where  n(»  is  commonly  referred  to  as  the  real  part  of 
the  refractive  index  and  k  is  the  absorption  index  of  the  contaminant  film.  Knowing  n(v)  and 
k|»,  one  can  calculate  the  transmittance  or  reflectance  of  any  contaminant  film  substrate 
for  a  given  film  thickness  through  the  use  of  Fresnel’s  equations  within  a  mathematical 
model.  Conversely,  n  and  k  can  be  theoretically  calculated  from  a  film-substrate 
combination  provided  transmittance  data  for  two  film  thicknesses  can  be  obtained  and  used 
within  a  mathematical  model  such  that  there  are  two  equations  involving  n  and  k,  the  two 
unknowns.  This  was  the  technique  utilized  in  determining  n  and  k  from  transmittance 
measurements  in  Refs.  1  through  4. 

In  measurements  of  contamination  of  any  facility,  a  method  for  determining  effects  of 
contaminant  as  a  function  of  contaminant  thickness  is  required.  Since  most  problems  are 
associated  with  films  a  few  microns  thick  or  less,  conventional  techniques  for  thickness 
measurement  do  not  apply.  However,  thin  film  optical  interference  provides  an  excellent 
technique  for  determining  thin  film  thickness,  e.g.,  film  thickness  from  0.25  to  25  /tm.  This 
requires  knowledge  of  the  refractive  index  at  some  known  wavelength.  After  having  found 
the  refractive  index  and  thickness,  there  are  techniques  available  for  determining  the 
complex  refractive  indices  (n’s  and  k’s)  of  the  film.  The  subject  of  this  report  is  the 
development  of  the  mathematical  models  and  computer  programs  required  for 
determination  of  n  and  k  from  transmittance  and  reflectance  experimental  data.  Computer 
programs  TRNLIN  and  SKKTRANS  are  used  for  analyzing  transmittance  data  of  multiple 
films  and  a  single  film,  respectively,  whereas  RENLIN  and  SKKRFL  utilize  reflectance  data 
of  multiple  films  and  a  single  film,  respectively. 

1.1  ANALYSIS  OF  TRANSMITTANCE  MEASUREMENTS  OF  THIN  FILMS 

When  a  beam  of  monochromatic  light  is  incident  upon  a  thin  film  deposited  on  a  thick 
substrate,  a  reduction  in  the  intensities  of  the  transmitted  and  reflected  beams  is  expected,  as 
well  as  phase  shifts  of  the  latter  two  beams  from  that  of  the  incident  beam.  At  a  given 
wavenumber,  v,  the  reduction  of  beam  intensities  and  the  phase  shifts  depend,  in  part,  upon 
the  values  of  the  film  refractive  index,  n(v),  and  its  absorption  index,  k(r). 
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In  one  method  of  obtaining  the  complex  refractive  index  ti(*»)  =  n(v)  +  ik(y)  of  dielectric 
films,  the  differences  that  changes  in  film  thickness  make  upon  the  absolute  transmittance 
are  observed  (Refs.  1  through  4).  (The  circumflex  denotes  a  complex  quantity.)  Because  only 
the  transmittance  is  observed,  at  least  two  separate  pairs  of  transmittance  and  thickness 
measurements  are  necessary  at  each  wavenumber,  v,  because  two  optical  parameters  of  the 
film  are  sought,  n(*>}  and  k(*0. 

If,  however,  the  phase  shifts  corresponding  to  each  transmittance  value  are  known,  then 
the  transmittance  measurement  for  only  one  value  of  the  film  thickness  is  needed  to  find  n(v) 
and  k(p).  Also,  the  assumption  upon  which  the  method  is  based,  that  the  optical  properties 
of  a  dielectric  film  are  independent  of  film  thickness,  is  not  required. 

Section  2.0  describes  a  method  for  finding  the  optical  constants  from  the  transmittance 
spectrum  of  a  single  thin  film,  even  when  one  cannot  directly  measure  the  change  of  phase 
spectrum.  In  it  is  discussed  the  pertinent  theory  and  the  computer  algorithms  derived  from 
them.  Also  shown  are  the  results  of  a  preliminary  investigation  of  some  existing  AEDC 
transmittance  data  for  20-K  cryofilms. 

1.2  ANALYSIS  OF  REFLECTANCE  MEASUREMENTS  OF  THIN  FILMS 

In  this  report,  two  methods  are  described  for  using  reflectance  measurements  to  find  the 
optical  constants  of  a  thin  dielectric  film  having  parallel  plane  interfaces  with  a  vacuum  and 
a  thick  metal  substrate.  The  reflectance  is  that  of  an  infrared  (IR)  beam  of  radiation  that  is 
nearly  normally  incident  upon  the  film  in  the  vacuum. 

In  the  more  common  of  the  two  methods,  at  a  given  wavenumber,  the  variance  in  the 
reflectance  corresponding  to  changes  in  film  thickness  is  observed.  This  procedure  is  called 
the  "multiple-film  thickness  method."  If  the  reflectance  measurements  of  at  least  two 
thicknesses  of  the  same  film  material  are  available,  the  experimenter  can  compute  the 
complex  index  of  refraction,  n  =  n  +  ik,  from  a  least-squares  algorithm. 

At  the  heart  of  this  method  is  the  assumption  that  the  optical  constants  n  and  k  of  a  film 
material  are  independent  of  the  film  thickness,  at  least  over  the  experimental  range  of 
thickness  values.  This  assumption  does  not  hold  at  all  for  very  thin  metallic  films  (Ref.  5), 
and  the  slight  dependence  of  n  and  k  of  N2/CO2  films  on  20-K  germanium  (Ge)  upon  their 
thickness  (Section  2.3.3)  is  demonstrated. 

1 

The  second  method  finds  n  and  k  of  a  film  from  reflectance  measurements  of  a  single¬ 
film  thickness,  and  the  assumed  independence  of  n  or  k  upon  the  film  thickness  is  not 
necessary.  With  this  method,  one  calculates  the  optical  constants  from  a  dispersion  relation 
relating  the  modulus  of  the  complex  reflection  coefficient  with  its  phase. 
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This  “single-film  thickness  method”  has  worked  successfully  when  applied  to  the 
reflectance  measurements  of  a  dielectric  film  deposited  on  a  medium  whose  modulus  of  the 
complex  index  was  lower  than  that  of  the  film  (Ref.  6),  but  no  one  has  fully  tested  its 
applicability  to  other  experimental  situations.  When  analyzing  the  reflectance  of  a  dielectric 
film  placed  upon  a  medium  whose  modulus  of  the  complex  index  is  higher,  one  runs  into 
many  more  mathematical  and  computational  difficulties  (Ref.  7).  The  way  these  problems 
are  attacked  is  the  subject  of  Section  3.0. 

In  Section  4.0  the  multiple-film  thickness  method,  which  is  incorporated  in  the  computer 
program  RENLIN,  is  described. 

2.0  EXTRACTION  OF  THE  OPTICAL  CONSTANTS  OF  A  THIN  FILM  FROM 
TRANSMISSION  MEASUREMENTS  OF  A  SINGLE-FILM  THICKNESS 

When  determining  the  complex  refractive  index  n  =  n  +  ik  of  thin  films,  most 
investigators  have  used  the  observed  dependence  of  transmittance  upon  film  thickness.  (See, 
for  example,  Refs.  1  through  4.)  However,  Maeda  et  al.  (Ref.  8)  have  employed  a  Kramers- 
Kronig  (KK)  dispersion  relationship  which  allowed  a  calculation  of  the  n-  and  k-values  from 
the  transmittance  spectrum  of  a  single-film  thickness.  These  authors  used  measurements  of 
the  transmittance  relative  to  that  at  zero-film  thickness.  Unfortunately,  errors  most  likely 
associated  with  estimates  of  quantities  outside  of  the  wavenumber  domain  of  their  data 
caused  very  noticeable  distortions  in  their  computed  n-  and  k-spectra. 

Described  in  this  section  is  a  new  subtractive  Kramers-Kronig  (SKK)  algorithm  which 
gives  highly  reliable  n-  and  k-spectra  from  measurements  of  the  absolute  transmittance  of  a 
single,  uniform  thin  film  deposited  on  a  uniform  thick  substrate.  Not  only  can  this  program 
reduce  the  number  of  necessary  experimental  measurements,  it  is  also  capable  of  noting  any 
differences  in  the  optical  constants  of  the  film  material  due  to  changes  in  film  thickness. 

2.1  THEORY 

Several  authors  (Refs.  1  through  5  and  8  through  10)  have  discussed  the  transmittance  of 
a  light  beam  through  a  thin  film  on  a  thick  substrate  where  all  surfaces  are  planar,  parallel, 
and,  ideally,  infinite  in  extent  (Fig.  1).  For  normal  incidence  of  a  semi-infinite  light  beam, 
the  absolute  transmittance  T(v)  at  wavenumber  v  is  given  by  (Refs.  1,  5,  and  8) 

T(v)  =  t*(»')t(»')  (I) 

where  t  is  the  complex  normal  transmission  coefficient  of  the  film  (medium  1  in  Fig.  1)  and 
the  substrate  (medium  2)  bounded  by  semi-infinite  media  0  and  3.  The  asterisk  indicates  the 
complex  conjugate  of  a  complex  quantity. 
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Figure  1.  Geometry  depicting  analytical  model  for  a  thin  film  formed 
upon  a  thick  substrate. 
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The  transmission  coefficient  is 

t  =  |  t23  |  toi2  exp  (  — ja2d2)/[l  -  |  r23  |2  j  r210  | 2  exp  (-  4&2&2)] w  (2) 

where  the  complex  Fresnel  transmission  coefficient  tjm  and  the  reflection  coefficients  rim  at 
the  interface  of  the  incident  (j)  and  refracted  (m)  media  are 

tjm  =  2nj/ (nj  +  nm)  (3) 

and 

^jm  =  (nj  —  Hm)/ (nj  +  Dm)  (4) 

for  normal  incidence.  In  Eq.  (2),  toi2  represents  the  complex  transmission  coefficient  for  a 
beam  going  from  medium  0,  through  the  film  of  thickness  di,  and  into  medium  2,  the 
substrate,  whose  thickness  d2  is  much  larger  than  the  coherence  length  of  the  beam.  The 
complex  reflection  coefficient  r210  refers  to  a  beam  incident  upon  the  film  within  the 
substrate.  These  coefficients  may  be  written  as 

toi2  =  toit|2  exp(i7idi)/[l  +  roir^  exp(2i7jd|)]  (5) 

f2io  =  foi  +  r,0  exp(2i-y1d|)]/[l  +  r2ir10  exp(2i7idi)]  (6) 

with 

02  =  2irvk2  (7) 

and 

7i  =  2irpn]  (8) 

Each  coefficient  accounts  for  the  infinite  number  of  reflections  within  the  film  which  cause 
thin  film  interference  effects.  The  quantity  02  arises  from  the  absorption  of  photons  within 
the  substrate,  whereas  7]  describes  the  phase  changes  as  well  as  the  photon  absorption  from 
beams  traversing  the  film.  Note  that  the  denominators  of  t<j,2  and  r2io  are  equal  because  rjTO 
=  -rmj  by  Eq.  (4). 

Equations  (1)  through  (8)  show  the  dependence  of  the  transmission  coefficient  t,  and 
thus  the  transmittance  T,  upon  eleven  parameters:  wavenumber  v ,  film  and  substrate 
thickness  d|  and  d2,  and  the  real  and  imaginary  parts  of  ho,  hi,  n2,  and  h3.  One  method  of 
finding  any  unknown  parameters  is  to  vary  at  least  one  other  parameter  in  a  known  manner 
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so  that  the  number  of  separate  measurements  of  T  at  a  given  wavenumber  is  at  least  as  great 
as  the  number  of  unknown  parameters.  For  example,  if  only  the  film  optical  constants  ni 
and  k]  are  unknown,  the  transmittance  T  must  be  measured  at  least  twice  at  every 
wavenumber  using  different  known  values  of  the  film  thickness  dL  (Refs.  1  through  4). 
Because  closed-form  solutions  for  ni  and  k|  from  Eqs.  (1)  through  (8)  are  not  possible,  the 
authors  of  Ref.  1  utilized  the  nonlinear  least-squares  algorithm  of  Marquardt  (Ref.  1 1)  in 
computer  program  TRNL1N. 

Alternatively,  the  minimum  number  of  experiments  may  be  reduced  to  half  that  required 
for  the  previous  methodif  it  is  possible  to  measure  t(v)  instead  of  just  T(p)  at  the  desired 
wavenumbers.  In  this  report,  values  of  the  optical  constants  nt  and  k(  of  the  film  are 
unknown  and  the  nine  other  parameters  in  Eqs.  (1)  through  (8)  are  known.  By  finding  t  at 
the  wavenumbers  of  interest,  experiments  may  be  done  on  just  one  Film  thickness  instead  of 
making  observations  of  at  least  two  different  film  thicknesses  when  only  T  is  obtained. 

The  equation 


t(v)  =  [T(p)]1*  exp[i0(p)]  (9) 

can  be  used  provided  that  <j>  (p),  the  change  in  phase  of  the  transmitted  light  beam  relative  to 
that  of  the  incident  beam  at  the  0-1  interface,  can  be  found.  [Note  that  Eq.  (2)  implies  that  t 
and  to  12  share  the  same  arcus,  or  phase  angle.]  Direct  measure  of  0  is  difficult,  expeciatly  for 
transmittance.  However,  far-infrared  asymmetric  Fourier-transform  spectroscopy  can  yield 
phase  information  for  reflected  beams  (Refs.  12  through  15). 

Even  if  the  phase  0(p)  cannot  be  directly  observed,  it  may  yet  be  determined  from  any 
one  of  several  equivalent  relationships  including  the  widely-used  Kramers-Kronig  (KK) 
dispersion  relations,  the  Fourier-allied  integrals  (Refs.  16  and  17),  and  the  Fourier-conjugate 
series  (Ref.  18).  In  this  work,  a  modified  form  of  the  KK  dispersion  equations  known  as  the 
subtractive  Kramers-Kronig  (SKK)  relations  (Refs.  1  and  19)  is  used.  Ahrenkiel  (Ref.  19)  has 
shown  that  errors  associated  with  the  extrapolations  which  must  be  applied  to  the  data 
because  of  its  finite  extent  are  reduced  in  the  SKK  formulation. 


2.1.1  Determination  of  Phase  0  of  t 

Following  Maeda  et  al.  (Ref.  8),  the  model  described  in  Eqs.  (1)  through  (8)  and  pictured 
in  Fig.  1 ,  is  the  KK  relation 
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w  ^  f®  hi  |t(,')|  dv'  ,  -  . 

+  2jx  =  -  P  j  0  -j— 5  +  2x«ii 


ip  f°°  fajV) 
IT  0  __  „2 


dv'  +  2irvd\ 


(10) 


where  j  is  an  integer  and  P  indicates  that  the  Cauchy  principal  value  of  the  integral  is  to  be 
used.  The  value  of  the  integer  j  is  immaterial  if  t(»0  is  computed  by  Eq.  (9),  However,  its 
value  is  not  arbitrary  and  must  be  known  where  the  SKK  relations  are  used. 

2.1.2  Derivation  of  the  SKK  Relation 


If  at  v  =  vq  the  optical  constants  of  the  film  material  under  study  are  obtained  from  a 
prior  measurement,  then  measurements  of  the  optical  constants  at  other  wavenumbers  may 
be  based  upon  those  at  v  =  vQ  in  the  following  manner.  At  v  -  v0,  Eq.  (10)  is 


0(1*0)  +  2m  x  = 


,0  *  fit  T(y') 

x  1  0  „-2_„02 


dv'  4-  2xpodi 


(11) 


where  the  integer  m  is  not  necessarily  equal  to  j  in  Eq.  (10).  Combining  Eqs.  (10)  and  (11) 
gives 

»  ftiT(v') 

0(„<2_  p2)(p'2  _  ,o2) 


4>(v)  +  2j  x  <Kvo)  +  2mx 


2  2 
V(f  -  Xs 


‘’O 


p ! 


or 


ftiT(f') 


p  p(p(j2  _  j,2j  w 

<Kv)  +  2j X  =  —  [0(*o)  +  2m x]  +  - P  J  Q  ~7~ ,  2  2 M  -2  2\ 

VQ  X  0  (V  i  -  ?ov 


dv' 


To  make  <p(v)  a  function  of  v  (i.e.,  a  single-valued  relation),  j  must  equal  m  because  at  v 
=  vq,  the  last  equation  is  <Kv 0)  +  2jx  =  0(*>q)  +  2mx.  Thus, 

V  r  -I  K^O2  -  "2)  f"  hi  T(p‘) 

«»)  +  2m* - «»o)  +  2m*  +  - p  L-rz - jj-75 - jr*'" 

VQ  x  0  (p'z  -  -  vo/ 

(12) 


Usually  the  values  of  n  and  k  of  the  film  at  i>0  are  known  from  previous  experiments. 
Equations  (2)  through  (8)  are  used  to  compute  t(»*o)  and  then  find  <j>(v0)  from 
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</»(po)  =  tan- 1  {lm[t(^o)]/Re[t(i'o)] }  (13) 

where  tan-1  (corresponding  to  subroutine  DATAN2  in  FORTRAN  IV)  has  a  range  interval 
of  (  -  -k,  v)  and  selects  the  value  that  lies  in  the  same  quadrant  in  which  t(j>0)  lies  when 
drawn  in  the  complex  plane. 

2.1.3  Determination  of  the  Integer  m  in  Eq.  (12) 


In  order  to  use  the  SKK  relationship  in  Eq.  (12),  the  integer  m  must  be  evaluated,  which 
by  Eq.  (1 1)  is 


m 


*o  r  faT(y') 

2a2  ^0  y'2  -  y02 


dp' 


<M»o) 

2r 


+  Mi 


(14) 


While  the  numerical  evaluation  of  the  integral  in  Eq.  (14)  is  straightforward,  it  is  nonetheless 
tedious  and  lengthy  because  of  the  difficulties  involved  with  calculating  the  contributions  to 
the  integral  from  the  wavenumber  domain  near  the  singular  point  v’  =  pq.  Fortunately,  the 
integral  term  is  almost  always  much  smaller  than  the  algebraic  sum  of  the  other  two  terms  in 
Eq.  (14)  and  may  be  safely  approximated  by  the  integral  term. 

Let 

N 

to  T(P)  *  to  T0  -  .E  aj  b(p  -  p) 

in  which  To  represents  an  average  value  of  the  absolute  transmittance  over  those 
wavenumber  regions  far  away  from  an  absorption  band,  *N  is  the  number  of  significant 
absorption  bands,  vs  is  the  wavenumber  at  the  minimum  transmittance  of  the  j,h  band,  and  6 
is  the  Dirac  delta  “function.”  The  constants  aj  are  aj  =  to  |To/T(vj)  | ,  so  now 

to  T(v)  =  to  T0  +  .£  to  [T(vj)/T0]  5{p  -  v j)  (15) 


The  integral  of  Eq.  (14)  becomes 


fttT(»')  “  to  T0  .  ,  n  £  MTC^/ToM*'  -  »j>  .  , 

P  )  — - -  dp'  *  P  i  - rdv'  +  P  I  £  - — - ; - —dp' 

J  0  y'2  -  PO2  0  y'2  -  r02  0  j  =  l 


v'2  -  VQ1 


The  first  integral  on  the  right-hand  side  is  zero  because  to  T0  is  constant,  and 


pr„-7 


dp’ 


0  p'2  - 
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is  zero.  The  second  integral  becomes 


fn[T(vj)/T0] 


giving  as  an  approximation  to  Eq.  14, 

vQ  N  fli[To/T(pj)] 

m  «  — —  L  — - - — 

2-x2  J  =  1  ~ 


0("o)  + 
2r 


^od] 


(16) 


In  the  present  study,  T0  is  chosen  to  be  T(^o)- 

2.1.4  Estimation  of  the  Complex  Index  of  Refraction  n  of  the  Film 

Many  times  in  the  course  of  extracting  the  complex  index  of  refraction  n(j>)  =  n(p) 
+  ik(y)  of  the  film  from  the  value  of  the  complex  transmission  coefficient  t(v),  it  is  found 
that  several  different  pairs  of  n-  and  k-values  will  yield  identical  t(p)  values  when  each  pair  is 
substituted,  in  turn,  into  Eqs.  (2)  through  (8).  In  order  to  select  the  “best”  pair  of  values  on 
physical  grounds,  it  is  necessary  to  resort  to  a  cruder  model  of  the  experimental  setup  from 
which  approximate  values  n  and  k  are  obtained. 


Near  strong  film  absorptions,  beams  traversing  the  film  are  rapidly  attenuated.  As  a 
consequence,  nearly  all  of  the  photons  which  have  frequencies  that  are  readily  absorbed  by 
the  film,  but  nonetheless  appear  in  the  transmitted  beam,  have  passed  through  the  film  just 
once.  If  all  of  the  photons  in  the  transmitted  beam  have  traversed  the  film  only  once,  and  if 
I3/I0  is  the  observed  ratio  of  the  transmitted  beam  intensity  I3  to  the  incident  beam  intensity 
I0,  an  estimate  of  k  is  (Ref.  1) 

k  -  fh  [  |  hj/fio  |  Ts/(I3/Io)]/(4r^d1)  (17) 

where  the  transmittance  Ts  of  the  bare  substrate  is 

Ts  =  t*sts  (18) 

and 

t»  =  1 123 1  to2  exp(  -  /32d2)/[  1  -  I  r23 1 2 1  r20 1 2  exp(  -  4/3^]  *  (19) 

[Eq.  (19)  can  be  heuristically  obtained  from  Eq.  (2)  by  elimination  of  the  “1”  subscript 
denoting  the  film.]  Often  one  can  put  measured  values  of  Ts  into  the  relation  indicated  in 
Eq.  (17). 
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Estimates  of  n  are  based  on  an  assumed  film  absorption  bandshape  as  a  Cauchy 
(Lorentzian),  Gaussian,  Voight,  or  similar  profile.  Peiponen  (Ref.  20)  has  derived  an 
expression  for  n  because  of  an  absorption  band  having  the  Cauchy  form 

k(*)  =  kan/fl  +  (4/y2)(v  -  i/,)2]  (20) 


Here,  kmax  is  the  value  of  k  at  v  =  v\,  where  v\  is  the  wavenumber  at  the  maximum 
absorption  and  7  is  the  difference  between  the  wavenumbers  at  the  half-maximum  points, 
i.e.,  7  is  the  full-width-at-half-maximum  intensity  (Fig.  2). 


If  all  N-bands  in  a  spectrum  are  assumed  to  have  the  Cauchy  shape  of  Eq.  (20),  an 
approximate  expression  for  n  is 


n(*0 


it  kmaxj  ''^Lj  (**  ^(’'Lj  —  •’j) 

n°  +  1-1  -  »j)2  +  ^(•'Lj  -  •'j)2 


(21) 


Here,  n0  is  the  contribution  to  n  because  of  absorption  outside  the  wavenumber  domain  of 
the  data.  For  the  jth  band,  i»j  is  the  wavenumber  at  which  the  maximum  absorption,  kmax  j, 
occurs,  and 

*Lj  ~  v)  ~  7/2  (22) 


is  the  smaller  half-maximum  intensity  wavenumber,  where  7j  is  the  difference  between  the 
larger  half-maximum  intensity  wavenumber  and  V]_-y  In  this  study,  no  was  equated  with 
n(M0),  and  the  other  parameters  were  estimated  in  Eq.  (21)  from  the  approximate  k-spectrum 
created  by  the  application  of  Eq.  (17)  to  all  measured  T(p)  values. 


2.1.5  Extrapolations  of  the  Data 


A  major  problem  in  using  the  SKK  relation  for  in  Eq.  (12)  is  that  the  integral  over 
the  wavenumber  in  Eq.  (12)  has  the  limits  zero  and  infinity.  It  is  necessary  to  extrapolate  the 
data  to  these  end  points,  and  one  now  investigates  the  asymptotic  behavior  of  fti  T(r)  for  the 
two  limits  v  —  00  and  v  —  0. 


As  v  approaches  infinity,  the  electrons  in  any  material  become  more  nearly  “free,”  and 
the  dielectric  constant  may  be  written  (Refs.  21  and  22)  as  e(y)  =  n2  (?)  »  1  -  vp2/p2, 
which  means 

n(f)  *  -  vp2/(2v2)  (23) 

where 

vp  =  o?p/(2irc)  =  2tjeV(mec)  (24) 
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Figure  2.  Refractive  and  absorption  indices  of  an  absorption  having 
a  Cauchy  (Lorentzian)  bandshape. 
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a>p  is  the  “plasma  frequency”,  17  is  the  number  density  of  free  electrons,  m*  is  the  mass  of  a 
free  electron,  e  is  the  electronic  charge,  and  c  is  the  speed  of  light  in  vacua.  It  follows  from 
the  relations  given  in  Eqs.  (I)  through  (8)  that  as  v  approaches  infinity, 

tjmW  *  1  +  (»W  -  v^)/{Av2) 

rjm(y)  *  (j'pm2  -  *'pj2)/(4>'2)  =  tjm(v)  -  1 

yielding 

|t012|2  =  t012*  toi2  »[l  +  Pp22/(4^2)  +  .  •  .]2 

I  ?210 1 2  =  ^210*  ^lO  =*  [2»»pi2{>,pl2  “  "p22)  +  vp2*  +  ■  •  .]/(16*^) 


and,  finally, 

T  =  t*  t  *  1  —  J'pi2/*8*'4)  +  •  •  •  or,  fti  T  *  -  j'P24/(8>'4)  (25) 

Here,  vpi  is  the  critical  “plasma  wavenumber”  defined  in  Eq.  (24)  for  the  jth  medium,  and 
Ppo  =  rP3  =  0  because  in  this  experiment  these  media  are  essentially  vacuums. 

As  v  approaches  zero,  for  dielectrics  (Refs.  21  and  22), 


e(»0  =  n2(»0  *  eo 


(26) 


a  constant,  whereas  for  conductors 

e(y)  =  n2(„) 


2ic 

Qv 


(27) 


where  a  is  the  conductivity  for  constant  currents.  For  the  dielectric  media  of  Fig.  1  from  Eq. 
(26), 

to  T(p)  *  fn  T0  (28) 

becomes  a  constant,  as  v  approaches  zero. 

Equations  (25)  and  (28)  are  used  to  calculate  the  end  contribution,  A  and  B,  of  the 
integral  in  Eq.  (12).  If  is  the  smallest  wavenumber  in  the  data  domain,  and  assuming 


fn  T(»>)  *  fn  T(><f) 


(29) 
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the  contribution  A  to  the  integral  over  the  interval  (0,vf)  is  given  as 


A  = 


v(vo2  _  p2) 


pj7 


hi  T(f') 


0  (y‘ 2  -  J'o2)(»''2  “  "2) 


d*' 


17  +  v 
Vf  -  v 


±to|-!L!Ji. 

Po  'pi  -  VO 


I) 


(30) 


Likewise  if  is  the  greatest  wavenumber  in  the  data  domain,  assume,  following  Eq. 
(25), 


to  T(»0  =  W/v4)  to  T(^ 


(3D 


in  the  interval  (?„,«).  The  end  contribution  B  over  (j»u,oo)  is 


H>o2  -  v2)  p  r  “ 

IT  \  ( V '2 


enT(p') 

vq2)(v‘2 


v(vo2  ~  v2)  "u4  to  Tft-u) 


[PQ2  +  p2  + 
vu  p04  v4 


I 

3  pu3  i>o2  v2 


+ 


1 

t'04(*'02  -  P2) 


Pu  -  p0  J 


(32)  ‘ 


2.2  COMPUTER  PROGRAM  SKKTRANS 

The  FORTRAN  IV  computer  program  called  SKKTRANS  extracts  the  optical 
parameters  n(y)  and  k(p)  of  a  film  material  from  the  absolute  transmittance  data  of  a  single 
fixed  film  and  substrate  surrounded  by  a  vacuum.  The  experimental  geometry  was  assumed 
to  be  that  of  Fig.  1  with  the  physical  description  contained  in  Eqs.  (1)  through  (8).  It  was 
further  assumed  that  the  user  would  supply  a  pair  of  n-  and  k-values  which  had  been 
obtained  in  a  separate  experiment  at  v  =  v0,  a  wavenumber  that  is  not  at,  or  next  to,  the 
ends  of  the  wavenumber  domain  of  the  data.  In  this  program,  then,  no  =  =  1  by 

assumption,  and  the  user  must  enter  the  rest  of  the  parameters:  film  and  substrate 
thicknesses  d}  and  and  the  refractive  index  n2  of  the  substrate. 
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The  program  SKKTRANS  is  divided  into  three  main  segments,  “MAIN,” 
“KAMKON,”  and  “NLNSYS.”  In  MAIN,  the  computer  reads  the  title  card,  job 
parameters,  and  the  data,  which  are  either  on  cards  or  in  temporary  direct-access  storage  if 
the  data  were  originally  on  tape.  (The  program  reads  tape  data  and  places  it  into  storage  in 
the  jobsteps  prior  to  MAIN.)  The  data  need  not  be  equally  spaced  in  wavenumber. 

The  second  and  third  parts  of  SKKTRANS  calculate  the  optical  constants  of  the  film, 
one  wavenumber  at  a  time.  KAMKON  calculates  t(v)  and  passes  it  on  to  NLNSYS  which 
finds  n(>0  and  k(p). 

2.2.1  Program  Segment  KAMKON 

In  KAMKON,  the  computer  recovers  the  necessary  phase  information  by  computing  the 
change  of  phase  <j>{v)  of  the  transmitted  beam  using  the  SKK  relation  [Eq.  (12)]  and  the 
integer  value  for  m  that  is  closest  to  the  right-hand  side  of  Eq.  (16).  The  computer  then  finds 
the  value  of  the  complex  transmission  coefficient  t(v)  using  Eq.  (9). 

The  computer  algorithm  in  KAMKON  which  evaluates  the  SKK  integral  of  Eq.  (12)  is 
very  similar  to  algorithms  in  other  SKK  programs  which  previously  were  supplied  to  AEDC 
(Ref.  1).  KAMKON  performs  the  integration  by  using  Simpson’s  rule  except  near  the 
singular  points  v'  -  vq  and  v'  =  v,  or  whenever  the  wavenumber  separations  between  three 
successive  data  points  are  unequal.  In  these  exceptional  cases,  the  numerator  of  the 
integrand,  fti  T(p),  is  approximated  by  a  quadratic  function  of  v  fitted  to  the  data  point  at 
the  wavenumber  of  interest  and  the  two  surrounding  data  points.  The  computer  uses  an 
analytical  expression  to  evaluate  the  contribution  to  the  integral  for  these  three  data  points 
and  computes  the  end  corrections  from  Eqs.  (30)  and  (32). 

2.2.2  Program  Segment  NLNSYS 

NLNSYS  uses  the  value  of  t(p)  that  is  found  in  KAMKON  and  calculates  the  film  optical 
constants  n(v)  and  k(v).  It  is  not  possible  to  solve  for  the  complex  refractive  index  m  in  Eqs. 
(2)  through  (8),  and,  consequently,  NLNSYS  incorporates  a  nonlinear  iterative  interpolation 
scheme  very  similar  to  Marquardt’s  nonlinear  least-squares  algorithm  (Ref.  11).  However, 
the  converged  values  of  the  optical  constants  which  the  present  algorithm  provides  cannot 
meet  the  least-squares  criterion  because  the  algorithm  employs  two  known  quantities,  the 
real  and  imaginary  parts  of  t(j»),  to  find  two  unknown  quantities,  n(j<)  and  k(»<). 

2.2.2.1  Selection  of  the  “Best”  Pair  of  Optical  Constants 

Always  from  strong  absorptions,  the  algorithm  in  NLNSYS  almost  always  converges  to 
the  “best”  pair  of  n-  and  k-values,  when  the  initial  guesses  of  these  values  are  the 
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corresponding  “best”  values  from  the  calculation  just  prior  to  the  present  one.  In  spectral 
regions  far  from  strong  absorptions,  the  “best”  pair  of  values  is  the  pair  in  which  n  is  not 
too  far  from  unity  and  k  is  close  to  zero.  The  other  pairs  of  n-  and  k-values  which  satisfy 
Eqs.  (2)  through  (8)  are  usually  not  close  to  the  “best”  pair  of  values. 

If  the  absorption  is  strong  at  the  wavenumber  of  interest,  however,  one  may  not  be  able 
to  readily  distinguish  the  “best”  pair  of  n-  and  k-values  from  the  other  possibilities.  In  this 
case,  NLNSYS  uses  at  least  eleven  different  initial  guesses  for  n  and  k,  and  usually  finds 
more  than  one  distinct  pair  of  converged  n-  and  k-values.  NLNSYS  uses  the  approximate 
values  of  n  and  k,  given  in  Eqs.  (21)  and  (17),  respectively,  to  help  in  the  selection  of  the 
“best”  pair. 

Theoretically,  if  the  absorption  band  has  the  Cauchy,  or  Lorentzian,  shape  of  Eq.  (20), 
the  n-values  reach  a  maximum  very  near  the  wavenumber  v  -  at  the  lower  half-maximum 
point,  and  a  minimum  very  near  the  wavenumber  v  +  corresponding  to  the  upper  half¬ 
maximum  intensity  point  (Refs.  20  and  23).  Thus,  in  the  wavenumber  interval  (v- ,  ?+)  n 
should  increase  monotonically  with  increasing  v  and  vice  versa  (Fig.  2). 

Experimentally,  measurement  errors,  non-Cauchy  bandshapes,  and  the  use  of 
approximate  n-  and  k-values  often  cause  shifts  in  the  location  of  the  extremal  values  of  n. 
Therefore,  expect  to  find  the  experimental  maximum  n-value  within  the  interval  (v- ~6, 
v-  +  5)  and  the  minimum  within  the  interval  (i>+-5,  v+  +  5),  where  6  is  typically  3  cm-1  or 
so.  Within  the  two  intervals,  the  n-values  should  not  be  rapidly  changing,  and  NLNSYS 
.  chooses  as  “best”  the  pair  of  n-  and  k-values  closest  (in  the  n,  k-plane)  to  the  n-  and  k-values 
found  in  the  last  calculation. 

Outside  of  the  two  wavenumber  intervals  mentioned  above,  NLNSYS  usually  picks  as 
the  “best”  values  the  pair  closest  to  the  approximate  n-  and  k-values  for  that  wavenumber. 
The  only  exceptions  occur  when  the  wavenumber  of  the  calculation  is  within  the  interval 
(p_  v  +  and  the  n-value  chosen  as  “best”  by  the  above  criteria  is  less  than  0.9  of  the 
n-value  found  at  the  larger  wavenumber  of  the  last  calculation.  (In  SKKTRANS,  the 
calculations  proceed  in  descending  wavenumber.)  Because  in  this  interval  the  n-values 
should  increase  monotonically  with  decreasing  wavenumber,  NLNSYS  chooses  as  “best” 
the  pair  having  the  smallest  n-value  that  equals  or  exceeds  0.9  of  the  best  n-value  from  the 
previous  ‘calculation. 

2.3  PRELIMINARY  RESULTS  FROM  SKKTRANS 

At  the  time  of  this  writing,  some  minor  debugging  of  SKKTRANS  may  still  be  necessary. 
However,  tests  have  been  made  on  the  reliability  of  the  SKK  and  interpolation  algorithms 
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using  a  theoretical  transmittance  spectrum.  The  effects  of  varying  the  film  thickness,  using 
some  existing  AEDC  cryofilm  data,  also  have  been  noted. 

2.3.1  Test  of  the  SKK  and  the  Interpolation  Algorithms 

A  theoretical  CO2  cryofilm  transmittance  spectrum  was  used  to  test  the  algorithms  in 
SKKTRANS.  This  transmittance  spectrum,  extending  from  700  to  3700  cm-1,  was 
computed  by  substituting  the  n-  and  k-values  of  pure  solid  C02  cryofilms  at  20  K  (Ref.  24) 
into  Eqs.  (1)  through  (8).  This  spectrum  was  then  processed  by  SKKTRANS  in  an  effort  to 
recover  the  original  n-  and  k-values. 

Note  that  the  SKKTRANS  did  not  simply  reverse  the  calculation  process  involved  in 
computing  the  test  transmittance  spectrum.  This  can  be  done  only  when  t  is  known  instead 
of  T.  SKKTRANS  does  compute  t,  by  first  finding  <t>,  and  then  it  calculates  the  optical 
constants  in  an  iterative  scheme  that  is  essentially  a  trial-and-error  method.  Also  note  that 
the  computation  of  n  and  k  at  a  single  wavenumber  requires  the  transmittance  values  at  all 
wavenumbers  included  in  the  data  because  of  the  SKK  relation  for  <j>  [Eq.  (12)]. 

In  essence,  the  validity  of  the  numerical  approximation  of  the  SKK  relations  [Eqs.  (12), 
(13),  (30),  and  (32)]  was  tested.  Figure  3  shows  the  excellent  agreement  of  the  SKKTRANS 
n- values  with  the  original  n-values  even  near  the  intense  vj  band  of  12Co2-  Here,  the  complex 
refractive  index  of  the  film  is  the  least-squares  value  n  =  1.2220  +  iO  at  the  wavenumber  v0 
=  2700  cm  -  *,  and  the  C02  film  thickness  is  taken  as  2.071  jun.  The  data  points  included  in 
the  SKK  algorithm  extend  from  710  to  3690  cm  - 1  and  are  spaced  every  10  cm  - 1 ,  except  near 
significant  absorptions.  For  this  film  and  all  others  discussed  in  this  report,  the  substrate  is 
20  K  germanium  4  mm  in  thickness  with  a  complex  index  of  refraction  n2  =  ng  +  iO  where, 
correcting  Eq.  (14)  of  Ref.  1, 

ng(?)  =  A  +  BL  +  CL2  +  Di>-2  +  Er4 
A  =  3.880 

B  =  3.91707  x  10”9  cm2 
C  =  1.63492  x  10-17  cm4 
D  =  -  600  cm2 
E.  —  5.3  X  10®  cm4 

and 

L  =  (v-2  -  2.8  x  10~ 8  cm2)-1  (33) 
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Figure  3.  Refractive  index  of  solid  COj  on  20-K  germanium. 
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Differences  between  the  SKKTRANS  n-spectrum  and  the  original  are  discernible  only  near 
the  relative  extrema  of  the  n-spectra.  The  exceptionally  close  agreement  of  the  SKKTRANS 
k-spectrum  with  the  original  is  evident  in  Fig.  4,  and  it  was  concluded  that  the  SKK  model 
and  the  approximations  can  faithfully  reproduce  the  original  spectra. 

Figure  3  also  includes  a  plot  of  the  approximate  n-spectrum  (NAPPROX)  computed 
from  Eq.  (21).  These  approximate  values  vary  widely  from  the  other  two  spectra  because  the 
?3  band  of  12Co2  is  highly  asymmetric  and  does  not  fit  the  single  Cauchy  (Lorentzian) 
bandshape  assumed  in  the  calculation  of  the  approximate  n-spectrum.  Nevertheless,  the 
criteria  discussed  in  Section  2.2.2. 1  for  choosing  the  “best”  n-  and  k-values  work  well  even 
in  this  case  of  a  poor  approximation  of  the  n-spectrum. 

On  the  'other  hand,  observe  in  Fig.  4  the  good  agreement  of  the  approximate  k-values, 
from  Eq.  (17),  with  the  other  k-spectra.  As  expected  from  the  assumptions  in  the  derivation 
of  Eq.  (17),  this  agreement  is  relatively  closer  as  the  k-values  increase.  The  results  are  in 
contrast  to  the  conclusion  of  Maeda  et  al.  (Ref.  8),  that  the  approximate,  or 
“experimental”,  k-values  of  Eq.  (17)  are  highly  distorted  from  the  “true”  KK  or  SKK 
values  contained  in  Eqs.  (10)  and  (12). 

2.3.2  Variation  of  +  2mx  with  Film  Thickness 

When  Eq.  (16)  was  derived,  it  was  assumed  the  Cauchy  principal  value  term  in  Eq.  (11) 
was  small  enough  to  make  <£(j>0)  +  2  mir  nearly  linearly  dependent  upon  the  film  thickness 
di.  [0(vo)  +  2  mir  is  the  phase  change  of  the  radiation  in  the  transmitted  beam  having 
wavenumber  vq.J  SKKTRANS  determines  from  Eq.  (13)  and  uses  Eq.  (16)  to  find  m. 
Thus,  a  plot  of  the  SKKTRANS  values  of  0(t'o)  +  2  mir  versus  film  thickness  is  an 
independent  check  of  Eq.  (1 1)  and  the  assumption  of  linearity. 

In  Fig.  5  are  seen  the  0(yo)  +  2  mir  values  found  by  SKKTRANS  for  15  different 
thicknesses,  ranging  from  0.2633  to  4.476  /*m,  of  films  formed  when  a  gaseous  mixture  of  N2 
(74.7  mole  percent)  and  C02  (25.3  mole  percent)  impinges  on  a  20-K  germanium  substrate 
(Ref.  1).  To  a  good  approximation  the  plot  is  linear  and  passes  through  the  origin,  and  the 
regression  line  (dashed  line)  through  the  data  has  a  slope  of  1.69967  ^m*1.  In  Eq.  (11),  the 
value  of  2  ttv^vq  =  2144  cm-1)  is  1.34711  /un-1  which  has  a  relative  difference  of  20.7 
percent  with  the  experimental  slope.  Because  the  authors  of  Ref.  1  found  the  film 
thicknesses  from  the  interference  pattern  of  He-Ne  laser  beams  reflected  from  the  film,  it  is 
possible  that  a  systematic  error  entered  the  thickness  measurements.  In  this  instance, 
increasing  each  of  the  measured  thickness  values  by  a  factor  of  1.69967  /un- 7 1.347 11  /un_1 
=  1.26171  would  cause  the  experimental  slope  to  be  identical  with  the  theoretical  slope,  2 
7rj>0.  However,  other  tests  (Section  2.3.3)  do  not  support  the  above  changes  in  the 
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Figure  4.  Absorption  index  of  solid  CO2  on  20-K  germanium. 
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Figure  5.  Plot  of  phase  shift  4>  =  <£q  +■  2mr  at  =  2144  cm  1  of  N2/CO2 
(74.7%/25.3%)  films  at  20-K  germanium. 
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thicknesses,  nor  does  the  experimental  error  estimated  by  the  authors  of  Ref.  1  allow  so 
great  a  discrepancy.  At  any  rate,  it  appears,  at  least  for  N2/CO2  films,  that  2  Tyodi  is  the 
dominant  term  of  Eq.  (11). 

2.3.3  Changes  in  the  n-  and  k-Spectra  as  a  Function  of  Film  Thickness  dj 

As  the  thickness  of  a  film  changes,  the  transmittance  spectrum  changes  in  appearance. 
However,  if  it  is  assumed  that  the  films  are  uniform  and  homogeneous  in  composition 
regardless  of  thickness,  two  films  of  the  same  material,  but  having  different  thicknesses, 
should  yield  the  same  n-  and  k-spectra.  This  assumption  is  the  basis  of  the  least- squares 
method  of  analyzing  transmittance  with  thickness  to  yield  the  film  optical  constants. 
However,  the  assumption  is  not  true  for  very  thin  metal  films  (Ref.  5).  Program 
SKKTRANS  allows  one  to  examine  each  film  thickness  separately. 

Figures  6  and  7  present  the  optical  constants  of  the  N2/CO2  films  discussed  in  Section 
2.3.2.  For  each  thickness,  the  data  points  were  spaced  2  cm  - 1  apart  and  extended  from  602 
to  3750  cm-1.  The  n(yo)-  and  k(vo>- values  required  in  the  SKK  algorithm  were  1.2610  and 
zero,  respectively,  and  were  measured  at  v0  =  2144  cm-1. 

The  dashed  line  in  Fig.  7  represents  the  values  of  k  found  in  Ref.  1  by  a  least-squares  fit 
of  transmittance  to  thickness.  The  authors  of  Ref.  1  computed  values  of  n  by  an  SKK 
analysis  of  the  least-squares  k-values,  and  the  dashed  line  in  Fig.  6  indicates  these  n-values. 
These  n-  and  k-spectra  are  nearly  identical  to  those  of  the  smallest  thickness  (0.2633  /rm) 
even  though  the  least-squares  analysis  used  the  transmittance  data  of  15  different  film 
thicknesses  ranging  from  0.2633  to  4.4761  /im.  This  is  because  the  transmittance  at  a  given 
wavenumber  decreases  approximately  exponentially  with  increasing  film  thickness,  and  the 
least-squares  criterion  is  biased  in  favor  of  the  larger  measurements.  For  the  N2/CO2  films, 
for  example,  the  measured  absolute  transmittance  at  2350  cm- 1  was  approximately  0.3  for  a 
film  0.2633-/im  thick,  whereas  the  1.3165-jun  film  had  a  value  of  about  0.01.  If  the  least- 
squares  procedure  had  achieved  nearly  the  same  absolute  error  (e.g.,  0.01)  between  the 
computed  and  measured  transmittance  values  of  these  films,  the  relative  error  would  be 
about  30  times  greater  for  the  thicker  film.  This  means  that  one  can  disregard  transmittance 
measurements  that  are  near  zero  when  employing  the  least-squares  method. 

There  are  also  reasons  to  discount  nearly-zero  transmittance  measurements  when  one 
makes  a  KK  or  SKK  analysis.  In  Figs.  6  and  7  only  the  SKKTRANS  n-  and  k-spectra  of  the 
smallest  five  thicknesses  are  shown  because  the.peak  absorptions  in  the  ten  larger  thicknesses 
are  so  intense  that,  because  of  instrumental  drift,  the  Fourier-transform  spectrometer 
measured  the  smallest  absolute  transmittance  to  be  very  slightly  negative.  In  cases  of  this 
type,  a  new  “baseline”  for  the  transmittance  spectrum  can  be  established  by  adding  a  very 
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Figure  6.  Refractive  index  of  solid  N2/CO2  mixture  (74.7%/25.3%) 
on  20-K  germanium  computed  at  several  thicknesses. 
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Figure  7.  Absorption  index  of  solid  Nj/COj  mixture  (74.7%/25.3%) 
on  20-K  germanium  computed  at  several  thicknesses. 
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small  constant  to  every  original  transmittance  value  so  that  the  new  transmittance  spectrum 
will  not  have  any  negative  (or  zero)  values. 

It  is  not  clear,  though,  what  the  additive  constant  should  be.  The  establishment  of  the 
criteria  for  choosing  the  constant  is  very  important  because  relatively  small  adjustments  of 
the  additive  constant  can  cause  variances  in  the  smallest  transmittances  of  several  hundred 
percent,  or  more,  and  can  create  very  noticeable  differences  in  the  n-spectra  computed  by 
SKKTRANS.  Perhaps  it  is  possible  to  find  a  suitable  additive  constant  from  “sum  rules” 
akin  to  those  developed  by  Smith  and  Manogue  (Ref.  25)  for  reflection  spectroscopy. 

Note  in  Figs.  6  and  7  that  both  the  n-  and  the  k-spectra  of  N2/CO2  vary  slightly  and  in  a 
nonrandom  manner  with  film  thickness.  The  same  effect  occurs  even  more  dramatically  in 
thin  films  of  85%  N2/15%  NH3,  91%  CO/9% H20,  and  93%  Ar/7%  H20  condensed  onto 
a  20-K  germanium  substrate  (Ref.  1).  Only  further  investigation  can  uncover  whether  this 
trend  is  a  real  physical  effect  of  these  films,  or  a  computational  artifact  due  to  a  systematic 
error  in,  e.g.,  the  determination  of  the  film  thickness,  or  inadequacies  in  the  model.  As  a  test 
of  the  effects  of  a  systematic  error  in  film  thickness,  the  optical  constants  of  the  N2/CO2 
films  were  recomputed  using  the  SKKTRANS  thickness  values  which  were  1.26171  times 
those  given  in  Figs.  6  and  7.  (In  Section  2.3.2,  it  was  found  that  1.26171  was  the  ratio  of  the 
observed-to-theoretical  slopes  of  the  plot  in  Fig.  5.)  Even  though  the  new  spectra  were  not 
identical  to  the  original,  the  new  spectra  showed,  to  the  same  degree,  the  original  variance  in 
value  with  film  thickness.  At  least  in  this  instance,  the  variation  of  n  and  k  does  not  seem  to 
depend  on  film  thickness  errors.  These  examples  illustrate  the  usefulness  of  obtaining  the 
optical  constants  of  individual  films,  especially  for  films  such  as  those  containing  hydrogen- 
bonded  molecules  whose  optical  properties  may  depend  on  film  thickness. 

2.4  USING  PROGRAM  SKKTRANS 

SKKTRANS  can  be  used  to  extract  the  optical  constants  of  a  uniform  thin  film  from 
absolute  transmittance  measurements  of  the  film  on  a  thick  substrate  surrounded  by  a 
vacuum.  Experimentally,  an  absolute  transmittance  measurement  is  the  ratio  of  the 
transmitted  beam  intensity,  at  a  given  wavenumber,  to  that  of  the  incident  beam. 

2.4.1  Program  Parameters 

The  present  version  of  SKKTRANS  calculates  the  refractive  index  of  the  20-K 
germanium  substrate,  n2  =  ng,  from  Eq.  (33).  One  must  enter  into  SKKTRANS  the  film  and 
substrate  thicknesses  and  the  optical  constants  of  the  film  at  v  =  v0.  a  wavenumber  that  is 
not  at  one  of  the  endpoints  of  the  transmittance  data. 
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The  program  also  requires  one  to  supply  the  wavenumbers  within  the  data  domain  at 
which  to  start  and  stop  the  calculations  and  to  specify  how  many  of  the  data  points  used  in 
the  SKK  algorithm  are  to  be  skipped  between  successive  determinations  of  the  film  optical 
constants.  The  initial  guesses  of  n  and  k  for  the  first  calculation  and  the  maximum  number 
of  iterations  in  the  interpolation  procedure  must  also  be  given. 

If  the  data  are  on  a  tape  or  disk,  SKKTRANS  constructs  a  new  data  set  from  the  original 
data.  The  smallest  and  largest  wavenumber  must  be  entered  in  the  new  data  set  and  the 
number  of  original  data  points  on  the  tape  or  disk  that  are  to  be  skipped  when  the  program 
picks  out  the  data  points  to  be  included  in  the  new  set  of  data.  Additional  program 
parameters  enable  SKKTRANS  to  write  the  computed  n-  and  k-values  onto  a  tape.  Please 
note  that  appropriate  changes  in  the  job  control  language  (JCL)  must  be  made  for  use  of 
tape  or  disk  volumes. 

The  remaining  input  parameters  have  the  information  needed  to  find,  in  the  approximate 
k-spectrum  generated  by  Eq.  (17),  the  wavenumbers  at  the  peak  and  the  two  half-intensity 
points  for  each  major  film  absorption  band.  At  a  wavenumber  near  a  large  absorption, 
SKKTRANS  makes  multiple  optical  constant  computations  and,  subsequently,  picks  the 
“best"  values  in  the  manner  described  in  Section  2.2.2. 1.  Empirically  it  was  found  that 
absorptions  whose  peak  k-values  are  a  few  tenths  or  less  can  be  ignored. 

2.4.2  Spectra  with  Overlapping  Bands 

In  the  case  of  two  or  more  large  overlapping  bands,  one  must  be  careful  in  the  estimates 
of  the  peak  wavenumber  and  of  the  full-width-at-half-maximum  intensity  of  each  band. 
This  is  necessary  because,  for  each  band,  the  search  for  the  wavenumbers  at  which  these 
points  occur,  in  the  approximate  k-spectrum,  is  within  the  wavenumber  domain  that  is 
centered  at  the  estimated  peak  absorption  wavenumber  and  extends  on  either  side  of  it  to  the 
estimated  full-width-at-half-maximum  intensity.  One  should  not  allow  the  search  for  these 
wavenumbers  for  a  band  to  extend  into  a  wavenumber  domain  where  absorptions 
attributable  to  other  bands  are  noticeable.  A  complete  listing  of  the  program  input 
parameters  and  their  formats  appears  at  the  beginning  of  the  large  MAIN  segment  of 
SKKTRANS. 

3.0  DETERMINATION  OF  THE  OPTICAL  CONSTANTS  OF  A  DIELECTRIC 
THIN  FILM  ON  A  METAL  SUBSTRATE  FROM  REFLECTANCE 
MEASUREMENTS  OF  A  SINGLE-FILM  THICKNESS 

In  the  past,  most  observers  have  used  reflectance  measurements  of  films  to  compute  the 
optical  constants  n  and  k  of  a  film  material  by  noting,  at  a  given  wavenumber  v,  how  the 
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reflectance  varies  with  film  thickness.  (See,  for  example.  Section  4.0  and  Ref.  26.)  This 
method  is  especially  useful  when  the  optical  properties  of  the  film  material  are  known  not  to 
change  with  film  thickness,  and  the  reflectance  data  are  available  for  at  least  two  values  of 
the  film  thickness. 

If,  however,  the  optical  constants  depend  on  film  thickness,  or  when  reflectance  data  for 
just  one  film  thickness  are  available,  one  must  find  the  constants  in  another  way.  In  Section 
2.0,  a  dispersion  analysis  of  the  transmittance  of  a  single  dielectric  thin  film  on  a  thick 
substrate  is  given  which  yielded  the  n-  and  k-spectra  of  the  dielectric  and  showed  that  the 
n-and  k-spectra  of  N2/CO2  films,  on  20-K  Ge,  vary  in  a  nonrandom  way  upon  the  film 
thickness.  At  the  core  of  that  analysis  was  the  calculation  of  the  phase  change  0  suffered  by 
the  transmitted  beam  using  a  dispersion  relation  between  0  and  the  modulus  of  the  complex 
transmission  coefficient.  The  phase  change  is  of  interest  for  its  own  sake,  and  the  recovery 
of  the  phase  change  from  intensity  data  constitutes  the  phase  retrieval  problem  (Refs.  27 
through  30). 

A  corresponding  analysis  of  reflectance  data  has  some  complications  not  found  in  the 
transmittance  analysis,  primarily  because  the  reflectance  can  be  zero  under  certain 
conditions,  whereas  in  principle  the  transmittance  can  never  be.  The  possibility  of  zero 
reflectances  introduces  additional  terms  into  the  usual  dispersion  relations  which  increase 
the  computational  effort  needed  to  use  the  single-film  thickness  method  (Refs.  7  and  31 
through  35). 

3.1  THEORY 

Figure  8  depicts  the  geometry  of  the  theoretical  model.  A  semi-infinite  beam  of 
monochromatic  radiation  is  normally  incident,  in  a  medium  of  index  n0  upon  a  thin  film, 
index  nj.  The  film  is  sandwiched  between  the  semi-infinite  media  having  indices  n0  and  n2, 
and  the  interfaces  are  plane-parallel  and  infinite  in  extent. 

In  a  typical  experiment,  the  incident  medium  is  very  nearly  a  vacuum  (no  =  1),  and  the 
medium  of  index  n2  is  a  metal  substrate.  One  is  able  to  consider  the  metal  substrate  as  semi- 
infinite  because  |  n2 1  is  much  greater  than  unity  for  IR  and  visible  radiation,  and,  hence, 
essentially  no  reflected  waves  exist  in  the  metal  substrate  if  it  is  sufficiently  thick. 

The  geometry  of  a  typical  experiment  usually  does  not  allow  the  detection  of  beams  that 
are  actually  normal  to  the  interfaces,  and  a  researcher  must  use  beams  that  have  nearly 
normal  incidence.  If  the  maximum  angle  of  incidence  of  a  ray  in  the  beam  is  less  than  10  deg 
or  so,  an  observer  frequently  does  not  have  to  distinguish  between  the  components  of  the 
beam  which  are  parallel  to  the  plane  of  incidence  and  those  which  are  perpendicular  to  it 
(Ref.  36). 
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3.1.1  Absolute  Reflectance  Rot2 

At  a  given  wavenumber,  the  ratio  of  the  intensity  of  the  reflected  beam  to  that  of  the 
incident  beam  is  the  experimental  absolute  reflectance  R0i2-  Theoretically,  the  model  in  Fig. 
8  gives  R012  as 

R012M  =  rStfWW*)  (34) 

where 

roi2  =  [foi  +  fi2  exp(2i7,di)]/[l  +  r0i  r!2  exp(2i7id,)]  (35) 

is  Eq.  (6),  if  the  “0"  and  “2”  subscripts  are  interchanged.  At  normal  incidence,  the  complex 
reflection  coefficient  ?jm  for  radiation  incident  in  medium  j  at  the  interface  with  medium  m 
is  given  in  Eq.  (4),  71  is  defined  in  Eq.  (8),  and  dj  is  the  film  thickness  (Section  2.0  and  Refs. 
5,  7,  37,  and  38).  Please  note  that  some  authors  define  ?jm  as  the  negative  of  the  right-hand 
side  of  Eq.  (4),  and  the  two  definitions  differ  by  a  phase  of  ±  x.  (For  examples,  see  Refs.  6, 
7,  and  25.) 

3.1.2  Relative  Reflectance  R 

In  many  experiments,  it  is  often  easier  to  measure  the  reflectance  of  an  optical  system 
relative  to  a  standard  reflectance  rather  than  its  absolute  reflectance.  The  reflectance  R  of 
the  vacuum-film-metal  system  relative  to  the  vacuum-metal  system  reflectance  is  defined  as 

D  _  !r  _  1r/Io  _  R012 

k. - -  -  (36) 

!rs  IrS/Io  R02 

Ir  and  Irs  are  the  reflected  beam  intensities  from  the  vacuum-metal  systems,  respectively,  lo 
is  the  incident  beam  intensity,  and  R02  is  the  absolute  reflectance  of  the  vacuum-metal 
system.  Note  that  R  can  be  greater  than  unity.  Values  of  the  complex  index  n2  of  the  metal 
can  come  from  measured  values  of  R02,  and  how  this  might  be  done  is  illustrated  in  the  next 
two  sections. 

3.1.3  Absolute  Reflectance  R02  of  the  Vacuum-Metal  Interface 

For  the  case  of  a  bare  metal  substrate,  the  theoretical  model  yields 

Ro2(»0  =  ro2(v)ro2(v)  (37) 

for  the  absolute  reflectance  at  the  vacuum-metal  interface,  and  the  corresponding  complex 
reflection  coefficient  is 
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T02  =  (no  -  n2)/(rlo  +  (38) 

In  an  experiment,  an  observer  often  fits  the  experimental  R02  values  to  a  function  of  the 
wavenumber,  such  as  a  polynominal.  For  instance,  Arnold  et  al.  (Ref.  26)  found 

Ro2(v)  =  0.939  -  (4.25  x  10"5  cm>  (39) 

to  be  a  good  fit  to  the  empirical  data  for  the  absolute  reflectance  of  a  vacuum-aluminum 
interface  in  the  infrared.  There  are,  however,  certain  theoretical  restrictions  on  the 
mathematical  form  of  Ro2(»0>  which  Eq.  (39)  does  not  meet,  and  these  are  taken  up  in 
Section  3.2.2. 

/ 

3.1.4  Evaluation  of  n2  from  R02 

A  researcher  must  know  the  complex  refractive  index  of  the  metal,  n2,  before  the  film 
index  hj  can  be  found  from  reflectance  data  because  it  is  necessary  to  compute  the  reflection 
coefficient  of  the  film-metal  interface,  ri2,  in  any  scheme  for  getting  rij.  [See  Eqs.  (4)  and 
(35).]  One  method  for  obtaining  the  metal  index  n2  is  through  a  Kramers-Kronig  (KK),  or  a 
subtractive  Kramers-Kronig  (SKK),  dispersion  analysis  of  the  observed  vacuum-metal 
reflectance  spectrum,  Ro2(v).  Two  computer  programs,  KKR  and  SKKR,  have  been  written 
which,  respectively,  perform  KK  and  SKK  analyses  upon  a  reflectance  spectrum  to  extract 
the  optical  constants  of  the  reflecting  medium. 

Another  way  to  calculate  the  complex  index  n2  of  the  metal  comes  by  solving  Eq.  (38)  for 
n2,  which  yields 

n2  =  n00  “  ro2)/U  +  r02)  =  n0  [l  -  R&  exp(i/302)]/[l  +  R$  exp(i/302)]  (40) 

where,  from  Eq.  (37),  tqi  is 

ro2  =  l  r02 1  expOftfl)  =  R&  exp(i/S02)  (41) 

The  phase  /S02  is  the  difference  between  the  phase  of  the  beam  reflected  at  a  vacuum-metal 
interface  to  that  of  the  incident  beam  at  the  interface  and  is  often  difficult  to  determine 
empirically. 

If  002  cannot  be  measured  directly,  then  good  estimates  of  it  can  frequently  be  made. 
Because  |  n2 1  is  much  greater  than  |  no  | ,  which  is  unity  for  a  true  vacuum,  the  complex 
reflection  coefficient  ro2  is  nearly  - 1  making  002  approximately  equal  to  it.  This  has  led 
many  authors,  Arnold  et  al.  (Ref.  26),  for  example,  to  simply  set  002 t0  *  throughout  their 
data  range. 
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j(3<)2  is  not  a  constant  function  of  the  wavenumber,  and,  furthermore,  if  002  is  set  to  t  and 
no  to  unity,  Eq.  (40)  gives 


n2  =  n2  +  ik2  =  (1  +  R£)/(l  -  R£)  +  i0  (42) 

Physically,  this  means  k2  of  the  metal  is  zero,  and  the  metal  medium  would  not  absorb  any 
of  the  radiation  impinging  upon  it!  In  fact,  a  metal  is  a  very  highly  absorbing  medium  over 
most  of  the  IR  and  visible  spectrum,  and,  for  metals  such  as  aluminum  (Refs.  39  through 
42),  k2  often  exceeds  n2.  At  8330  cm-1,  for  instance,  bulk  A1  has  the  optical  constants  n2 
=  0.78  and  k2  =  9.16  (Ref.  41). 


The  major  problem  in  using  Eq.  (42)  as  an  estimate  of  n2  is  that  sometimes  this  can 
introduce  sizeable  errors  in  the  value  of  the  reflection  coefficient  at  the  film- metal  interface 

f,2  =  (nj  -  n2)/(ni  +  112)=  |r12j  exp(i/312)  (43) 

with 

012  =  tan-1  {[(n  +  n2)(k  -  k2)  +  (n2  -  n)(k  +  k2)] 

/[(n 2  -  nj)  +  (k2  -  k^)]j  (44) 

and  nl  =  n  +  ik  for  the  film.  If  k2  is  set  to  zero,  so  that  0O2  is  tt  and  Eq.  (42)  gives  n2,  the 
phase  change  upon  reflection  at  the  film-metal  interface,  0i2,  becomes  a  ficticious  value 

0{2  =  tan- *[2n2k/(n2  -  nj  +  k2)]  (45) 

In  the  spectral  regions  where  the  modulus  of  the  film  index  of  refraction,  |  ni  j ,  is  much  less 
than  that  of  the  metal,  |n2|,  both  0i2  [Eq.  (44)]  and  0|'2  [Eq.  (45)]  are  approximately  t  in 
value,  and  we  can  safely  use  Eq.  (42).  This  situation  is  often  encountered  experimentally, 
and  some  authors  set  012  as  well  as  0O2  to  *  throughout  their  data  domain  (Ref.  26). 

However,  if  either  of  the  optical  constants  of  the  film  of  the  metal  become  comparable  in 
magnitude  to  one  of  the  metal  constants,  some  larger  errors  are  possible.  Suppose  k  and  k2 
are  approximately  equal  in  a  region  of  strong  absorption  by  the  film,  but  n  is  very  small 
compared  to  n2.  The  true  value  of  0|2  is,  from  Eq.  (44),  012  *  tan-i[2k/(-n2)],  which 
means  0i2  is  about  0.578  tt  for  the  case  k  =  k2  =  2n2.  In  contrast,  if  002  is  set  to  t,  and  thus 
has  a  fictitious  k2  value  of  zero,  0i2  from  Eq.  (45)  is  about  tan  - 1  [2n2k/(k2  -  n^)],  giving  a 
fictitious  value  of  about  0.295  t  when  k  =  2n2.  The  discrepancy  between  these  two  values  of 
0i2,  0.283  it,  is  large  enough  to  cause  a  noticeable  error  in  r!2  and  roiz,  if  Eq.  (42)  is  used  to 
compute  n2.  Fortunately,  instances  in  which  k  and  k2  are  comparable  rarely  come  up  except 
in  spectral  regions  where  the  film  has  strong  absorptions. 
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There  is  a  greater  chance  of  error  if  one  sets  /3]2  as  well  as  002  to  it.  In  the  above  example, 
for  instance,  the  true  value  of  /?  12  is  closer  to  Vi  v  than  t. 

3.1.5  Multiple-Film  Thickness  Method 

The  model  depicted  in  Fig.  8  and  presented  in  Eqs.  (34)  through  (44)  is  the  basis  of  the 
method  for  finding  the  film  optical  constants  at  a  given  wavenumber  by  noting  the  change  in 
the  relative  reflectance  R  as  the  Film  thickness  varies.  This  multiple-film  thickness  method  is 
valuable  for  films  which  do  not  change  their  optical  properties  with  thickness  and  is 
discussed  more  fully  in  Section  4.0. 

3.1.6  Single-Film  Thickness  Method 

In  cases  where  the  optical  constants  of  the  film  change  with  its  thickness,  or  where  there 
are  reflectance  data  for  a  film  material  at  only  one  thickness,  one  cannot  apply  the  multiple- 
film  thickness  method  of  Section  4.0,  because  one  seeks  n  and  k  at  a  given  wavenumber  and 
only  one  experimental  piece  of  information  exists,  namely,  R. 

A  second  piece  of  information  can  be  supplied  from  the  use  of  any  of  the  appropriate 
dispersion  relations  between  the  modulus  of  a  complex  quantity  and  its  phase  (Refs.  16 
through  18, 22,  and  43).  Toll  (Ref.  44)  has  shown  that  postulating  the  existence  of  dispersion 
relations  for  a  linear  system  is  equivalent  to  assuming  that  strict  causality  (“no  output  before 
input”)  holds  for  the  system. 

In  the  present  instance,  one  can  find  the  modulus  of  the  reflection  coefficient  roi2  of  the 
vacuum-film-metal  system  directly  from  experiment  by  using 

|r0i2|  =  R0ft  =  (RR«2)/l  (46) 

where 

roi2  =  1 1012 1  exp(i^)  (47) 

where  $  is  the  phase  change  upon  reflection  at  the  vacuum-film  interface  for  a  beam  incident 
in  the  vacuum.  The  phase  change  <f>  is  what  is  sought  from  a  dispersion  relation  so  that 
optical  constants  of  the  Film,  n  and  k,  can  be  extracted  by  inversion  of  Eq.  (35). 

Lupashko  et  al.  (Refs.  6  and  7)  have  derived  the  Kramers-Kronig  (KK)  dispersion 
relations  for  the  special  case  in  which  a  planar  thin  film  is  sandwiched  between  two 
transparent  semi-infinite  media.  The  following  derivation  of  the  KK  relations  of  the 
vacuum-Film-metal  system  is  along  the  lines  of  Lupashko  et  al.  (Ref.  7). 
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3.1.7  Derivation  of  the  Kramers-Kronig  (KK)  Dispersion  Relations  for  roii 

In  order  to  proceed,  one  must  find  and  examine  the  analytic  extension  of  the  theoretical 
expansion  for  roi2  (?)  into  the  upper  half  of  the  complex  wavenumber  (i>)  plane.  A  function 
f(z)  that  has  a  derivative  relative  to  the  complex  number  z  at  all  points  within  a  region  of  the 
z-plane  is  an  analytic  function  of  z  within  that  region  (Refs.  45  through  47).  The  complex 
wavenumber  is  denoted  by 

V  =  +  i?!  (48) 


3.1. 7.1  Extension  of  Optical  Quantities  into  the  Complex  p-Plane 

Landau  and  Lifshitz  (Refs.  22,  35,  and  43)  have  demonstrated  that  in  the  upper  half  of 
the  p-plane,  the  extension  of  index  hj  for  material  j  has  the  following  properties: 

1.  The  hj  are  analytic. 

2.  The  values  of  hj  on  either  side  of  the  positive  imaginary  axis  are  paired  in  the 
following  way:  hj(  —  p*)  =  h*(?). 

3.  As  v  approaches  infinity  in  any  direction  in  the  upper-half  plane,  hj(v) 
approaches  unity. 

4.  Along  the  positive  imaginary  axis,  where  hj  is  real  (Property  2),  the  values 
decrease  monotonically  from  hj  (iO),  which  is  greater  than  unity,  to  hj  (i»)  =  ] 
(Property  3). 

5.  The  hj  are  real  only  on  the  imaginary  axis. 

From  a  straightforward  application  of  Property  2  to  the  extensions  into  the  i>-plane  of 
Eqs.  (4),  (8),  and  (35),  respectively,  it  is  found  that  in  the  upper-half  plane  (Ref.  43) 

*jm(  ~  £*)  =  fjm(^) 

and 

roi2C  —  =  foi2(*0  (49) 

3.1.7.2  Poles  and  Zeros  of  r0n  (v) 

Dispersion  relations  can  come  from  the  real  and  imaginary  parts  of  the  closed  contour 
integral  4  S(v)dv/(p2  -  v2)  (Refs.  7  and  35),  which  is  zero  if  S(>>)  is  analytic  in  the  region 
contained  within  the  closed  path  C  lying  in  the  upper  half  of  the  p-plane,  i.e., 
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$  S(»-)d£/(i2  -  v2)  =  0  (50) 

c 

if  S  is  analytic  within  path  C.  It  must  be  insured  that  C  circumvents  any  poles  of  S  as  well  as 
i  —  ±  v  on  the  real  axis. 

A 

In  this  case,  S  is  chosen  along  the  real  axis  to  be 

S(»)  =  ftirOi2(>0  =  fti  Jr0i2(r)|  +  (51) 

because  interest  is  in  the  relationship  between  the  modulus  of  roi2>  which  can  be  measured, 
and  its  phase.  The  analytic  extension  of  S  into  the  upper-half  plane  is  formally 

S(v)  =  fn  r0i2(»')  =  fti  |  r0i2  (v)  |  +  i^t*)  (52) 

Equation  (52)  indicates  that  any  pole  or  zero  of  f0i2(»i)  is  also  a  pole  of  S(£).  However, 
there  are  no  poles  of  r0i2  because,  in  the  extension  of  Eq.  (35),  j  fQi(v)f |2(v)  exp[2i7i  (v)dil ; 
is  less  than  unity  everywhere  in  the  upper-half  plane  (Ref.  7). 


The  zeros  £j  of  r0i2  (P)  can  be  found  by  setting  the  extension  of  Eq.  (35)  to  zero.  With  the 
aid  of  the  extension  of  Eq.  (8),  the  location  of  the  zeros  at  v  =  are  for 


=  fti[-roi(pj)/ri2(rj)]/[4«ini(vj)di] 
=  { Ln [  —  r0 i(Pj)/ r i2( ^j)]  +  i  •  2pirj 

/  j^ihi^dj 


(53) 


The  Ln  function  has  the  principal  value  of  the  natural  logarithm  function,  fn,  and  p  is  any 
integer.  Equation  (53)  is  transcendental,  so  it  cannot  be  solved  for  directly.  The  relations 
are  useful,  though,  in  an  iterative  search  for  v}  (Section  3.2.3). 


According  to  Eq.  (49),  if 


=  *Rj  +  i  K|j 


(54) 


is  a  zero  of  roi2,  -  ?j*  =  -  VRj  +  ifjj  is  also  a  zero;  that  is,  the  zeros  of  ?oi2  are  reflected 
across  the  imaginary  axis. 
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3.1.7.3  Evaluation  of  £  S(?)d£/(P  -  v2)  =  0 

The  closed  path  of  integration,  C,  (Fig.  9)  is  basically  a  quarter  circle  of  radius  L  and  the 
radial  lines  at  each  end  of  the  arc  which  are  along  the  positive  real  and  imaginary  axes.  The 
path  avoids  the  pole  at  v  —  v  on  the  positive  real  axis  with  a  semicircle  of  radius  and  goes 
around  a  possible  pole  of  S  at  v  =  irI0  on  the  positive  imaginary  axis  with  a  semicircle  of 
radius  go-  There  are  cuts  and  7/ and  small  circles  of  radius  gj  about  the  Z  zeros  of  roi2, 
not  on  the  imaginary  axis,  which  are  also  poles  and  branch  points  of  S.  This  choice  of  path  C 
makes  S(i»)/(i'2  -  v2)  analytic  within  it.  Thus,  Eq.  (50)  may  be  used,  which  is 


Figure  9.  Path  C  of  contour  integration  $  S(t^)dP/(v2  -  v2). 

C 
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q  _  ye  SC^iOd^R  [  |L  S(pR)dyR  |  0  S[k  +  QCxp(ia)]d[eexp(ia)] 

0  vfa  -  v2  r+Q  —  v2  a=  ir  [i>  +  gexp(ia)]2  —  i>2 


£  f  \  S(7j)d 

j  =  i  ^Ti  (VRj  +  ip 


S(7j)d(ifi)  J  S(y/)d(ii']) 


j)2  _  j/2  7j  (pRj  +  i„[)2  _  v2 


,  °j~°  S[pj  +  ej  exp(iaj)]d[gjexp(igj)]  1 
“j  =  2t  [v]  +  cjexp(io!j)]2  -  v2  ) 


f9~x/2  S[Lexp(ifl)3d[Lexp(ifl)] 

®=0  [Lexp(i0)]2  -  v2 

j L  S(iy[)d(jgi)  +  j0  SQpQdQyi) 

^io+eo  n2  +  v2  fIO-po  "i2  +  y2 


J1  ~*n  S[it*io  +  poexp(iqf0)]d[go  exp(ia0)] 
a0=*/Z  [ii'IO  +  CoexpOao)]2  -  y2 


(55) 


If  no  zero  of  f0i2  appears  on  the  positive  imaginary  axis,  the  last  three  integrals  on  the  right- 
hand  side  of  Eq.  (55)  are  combined  into  a  single  term: 


J 


L  S(ipj)d(ipj) 

0  i'!2  +  v 2 


Now  proceed  to  simplify  Eq.  (55)  by  reducing  q  and  all  gj  to  zero  and  letting  L  approach 
infinity. 


The  sum  of  the  integrals  may  be  written  along  path  segments  y}  and  yf  as  q}  goes  to  zero 


as 


J  S(7j)d(in)  1,  S(7j)d(i?i) 

^  (^Rj  +  i^l)2  ~  »2  ^  (•'Rj  +  ii'i)2  ~  *2 


*-lj  ASjdj/i 

=  -l  I 


0  (PRj  +  i^l)2  -  v2 
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-  "ij  (y|j  -  s>2  -  v\  -  2ivRjt»j)dyi 
J  0  p?  +  2  (4j  +  f2)  *l  +  ("Rj  -  l’2)2 


-  tan-1 


yij 

"Rj  +  * 


[•'I  +  (PRj  -  v )2  ^  (yRj  +  y)2  1  1 
•'Ij  +  (j'Rj  +  v)1  (j'lj  ~  I')2  i 


(56) 


The  quantity  ASj  =  S(7j)  -  S(7j)  is  the  value  of  S  at  a  point  on  y-  minus  the  value  of  S  at  the 
corresponding  point  on  7 j.  In  general,  ASj  is  not  zero  because  of  the  excursion  of  path  C 
about  the  branch  point  t>}t  although,  in  our  case,  ASj  is  a  constant  value  (see  Section  3. 1.7. 4) 
and  we  can  remove  it  from  the  integrand. 


In  the  limit  as  all  Qj  go  to  zero,  the  integrals  about  all  of  the  branch  points  of  S  are  zero, 
i.e.. 


f  S[pj  +  eiexp(iaj)]d[ejexp(io!j)]  ,  ^ 

lim  )  - = - -—7j - 3 -  =  0,  j  =  0,  1,  2,...,  Z 

cj  0  aj  [pj  +  ejexp(iaj)]2  -  p2 


(57) 


because  S  is  only  logarithmically  infinite  at  v  -  v\  (Ref.  35). 

The  limit  about  the  pole  v  =  v  is  taken  as  e  approaches  zero  and  L  gets  infinitely  large: 


r  S[p  +  eexp(ia)]d[eexp(ia)]  .  f0  S[p  +  eexp(i«)l 

H  J  - — - — - 5 - =  1  j  l«n_ - —  da 


[p  +  gexp(ia)]2  —  p2 


*■  c-o  2p  +  gexp(ia) 


i7rS(p) 

Iv 


(58) 


and  that  of  the  sum  of  the  first  two  integrals  in  Eq.  (55): 


lim 
e-o 
L— 00 


[/ 


9  S(pR)dPR  +  |  L  S(PR)dPR 


pi  -  p2 


p  +  e  pR 


2  -  p2 


1 


p ! 

0 


00  S(pr)  dpR 


(59) 


P  indicates  the  Cauchy  principal  value  of  the  integral. 
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Equations  (56)  through  (59)  are  incorporated  when  Eq.  (55)  is  rewritten  for 
infinitesimally  small  q  and  gj  values  and  infinitely  large  L  values  as 


-  r”  SQgjdTR 

o  4  -  "2 


iTrS(j') 

2v 


i  &  - 

—  ASj 
2v  j  =  i  J 


|^tan_1 


"ij 


PRj  -  * 


-  tan-1 


"ij 


‘'Rj  +  V 


—  ftl 
2 


r  vi j  +  (*,Rj  -  v)1 

(*Tlj  +  v)2  -I  1 

*-  »»Ij2  +  ("Rj  +  v)2 

(vRj  -  v)2  J  ) 

..  .  'n/1  S[Lexp(ifl)]dt?  .  t“  S(i*[)dj/, 

+  lim  i  ]  - +  i  )  - — 

L-0°  0  Lexpl(ifl)  -  p2/(Lexp(ifi)] )  0  P|2  +  v1  (60) 

3. 1.7.4  Case  S(^)  =  fti  ronM 

In  this  section,  Eq.  (60)  is  evaluated  for  the  choice  of  S,  namely, 

S(*)  =  fit  r0j2  (*) 

The  change  in  the  value  of  S  along  the  circular  paths  about  the  branch  points  Pj  is 

ASj  =  S(7j)  -  S(7j)  =  -  2iri  (61) 

for  j  =  1,  2,  ...»  Z.  The  negative  sign  indicates  the  clockwise  route  of  the  integration  path 
about  ij. 

The  integral  about  the  quarter  circle  of  radius  L  depends  on  the  limit  of  the  integrand  as 
L  goes  to  infinity.  Along  the  real  axis  the  asymptotic  behavior  of  index  n  of  any  material  is, 
by  Eq.  (23)  for  very  large  j»r, 

n(j'R)  *  -  +  *«*  (62) 

where  vp  is  the  “plasma  wavenumber”  defined  in  Eq.  (24).  Substituting  Eq.  (62)  into  Eqs. 
(4)  and  (8)  and  placing  that  result  into  Eq.  (35),  the  asymptotic  expression  of  ?0i2  for  large  pr 
is 

ioi2  (pr)  “  [^p?  +  (('p \  -  J'pf)exp(i4iTPRd1)  +  ...]/{44)  (63) 
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Here,  as  elsewhere  in  this  report,  the  subscripts  0,  1,  and  2  refer  to  the  corresponding  media 
in  Fig.  8. 

The  analytic  extension  of  Eq.  (63)  is  used  to  calculate  the  limit  of  the  integrand 
S[  Lexp(i  0)  ]  /  { Lex  p(  i0)  -  v2/[Lexp(i0)]  j  for  infinitely  large  L.  This  is  (Ref.  44) 

..  fa  rou[Lexp(i0)]  fatvpf/^L2)]  -  2i0 

lim  -  =  lim  - - - 

L“"  Lexp(i0)  L~ao  Lexp(i0) 

=  2exp(-  i0)  lim  L_1  fh  L_1 


=  0 

Hence,  the  integral  along  the  arc  at  infinity  is  zero,  i.e., 


t/2 


i 


S[Lexp(i0)]d0 
Lexp(i0)  -  p2/[Lexp(i0)] 


=  0 


(64) 


All  that  remains  to  evaluate  in  Eq.  (6)  is  the  integral  along  the  positive  imaginary  axis.  To 
do  this  one  must  know  the  behavior  of  roi2  on  the  imaginary  axis. 


From  the  discussion  in  Section  3. 1.7.1  it  is  known  that  n0,  n(v  n2,  r0i,  ri2>  and  y\  are  real 
everywhere  on  the  positive  imaginary  axis,  and  hence,  by  Eqs.  (8)  and  (35),  so  is  r0i2.  Also 
known  is  that  no  is  unity  everywhere,  and  ni(iv,)  and  n2(ivi)  are  greater  than  unity  for  finite 
positive  »i.  This  means 

- 1  <  roi(ij'i)  <  0,  if  0  s  vx  <  w  (65) 


On  the  other  hand,  r^O^i)  can  have  either  sign  depending  on  the  difference  rii(i^i)  —  h2(ivj). 


It  is  also  true  that 


exp[2i-y1(i»'i)d|]  =  exp[-  4Td1n,(iv1)^I]  £  1, 


(66) 


if  >  0. 

These  relations  are  sufficient  to  indicate  the  behavior  of  r0i2  on  the  positive  imaginary 
axis.  At  the  origin  (P  =  iO),  the  metal  index  h2  is  infinite  (Refs.  22  and  43),  and  r12(i0).and 
?oi2(iO)  are  both  equal  to  -  1,  i.e.. 
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?i200)  =  roi2(iO)  =  -  1 


(67) 


At  infinitely  large  vu  the  indices  approach  unity  (Property  3  of  Section  3. 1.7.1),  and 

*0120  00 )  =  0  (68) 


If  ni(i^i)  is  always  less  than  n^ivi),  and  r^iv)  and  roi2(ii»i)  never  become  positive,  then 
no  finite  zero  of  r0i2  can  exist  on  the  positive  imaginary  axis.  This  case  of  no  zeros  seems 
very  likely  because  while  ht  has  a  finite  value  at  the  origin,  n2  is  infinite  there,  and  both 
indices  decrease  monotonically  to  unity  as  approaches  infinity  (Property  4).  Even  if  hi 
should  become  greater  than  n2  somewhere  on  the  imaginary  axis,  making  r12  positive  at  that 
point,  r0|2  will  remain  negative  as  the  distance  is  increased  from  the  origin  until  r0j  is  the 
negative  or  ri2exp(-47ri'|dihi),  which  would  occur  at  the  first  zero  point,  v  =  ivi0.  Other 
zero  points  may  appear  even  further  from  the  origin  with  r0i2  alternating  in  sign  in  adjacent 
intervals  between  zeros.  In  any  case,  the  number  of  zeros  on  the  imaginary  axis  should  be 
very  few  because  n2  is  usually  greater  than  hi,  and  for  sufficiently  large  vu  the  exponential 
decay  will  dominate  r12  exp{  -  47rv(dih|). 


If  the  only  finite  zero  of  r0L2  on  the  positive  imaginary  axis  occurs  at  v  =  ivyo,  it  can  be 
written 


roi2(ivi)  =  ~  I*oi20*’i)l*  f°r  0  ^  Vy  <  P[0 


lroi2(i*'i)|»  for  vi  >  *[io 


(69) 


and  the  integral  in  Eq.  (60)  along  the  imaginary  axis  is 

.  P  r°°  _  p  f00  fri  |  r0i2(iyi)  |  di^i  _  ,  piq 

0  v\  +  v2  0  vy2  +  v2  v  v  (70) 

It  turns  out  that  if  the  only  zero  of  r0i2  on  the  imaginary  axis  is  the  one  always  present  at 
infinity  [Eq.  (68)],  then  the  same  answer  will  be  derived  as  when  j»[0  =  oo  is  substituted  into 
Eq.  (70).  In  that  instance,  the  arctangent  term  is  -  tt2/(2v).  If  there  are  two  or  more  finite 
zeros  of  foi2  on  the  positive  imaginary  axis,  each  zero  contributes  an  analogous  arctangent 
term  which  alternates  in  sign  in  ascending  v\. 

Now  the  desired  relations,  Eqs.  (61),  (64),  and  (70),  can  be  written  for  Eq.  (60)  for  the 
case  S  =  fnroi2.  When  these  relations  are  used,  Eq.  (60)  is 
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q  _  p  f  falroi2(xR)|dxR  +  r® 
0  *  ’■  1  - 


2  2 
VR  -  * 


0  4  -  *'2 


iirfti|r0i2(p)j  *tf(p) 


2p 


2p 


-  .E  [tan 

»  J-l 


-1 


"Ij 


XRj 


—  V 


-  tan 


-l 


•'ij 


XRj  + 


-  ftl 

2 


r  vij  +  (xrj  -  x)2 

(PR  +  p)21  1 

L  Vlj  +  (pRj  4-  v)2 

*  (*R  -  x)ij  ) 

+  IP  -  iun-1  J» 

0  v\  +  j>2  V  v 


(71) 


3.1.7.5  Kramers-Kronlg  Dispersion  Relation  for  <Kv) 


The  solution  of  the  real  part  of  Eq.  (71)  for  the  phase  #  of  f0i2  is  the  Kramers-Kronig 
(KK)  relation 

,  ~2y  n  7  ft,l?012(xR)|dPR 

0(x)  +  2m  ir  =  - P  j  - = - = - 

v  0  PR  -  v 


N  / 
.E  (t 
i=i  v 


+  2  E  ( tan-1 


•'ij 


XRj  -  V 


-  tan-1 


XRj 


_ \ 

Lj  +  V  ) 


-  2  tan  ~ 1 


pio 


(72) 


m  is  an  integer,  and  the  addition  of  a  multiple  of  2tt  to  Eq.  (72)  does  not  affect  exp(i</>). 
Zeros  of  r0i2  occur  at  v j  =  pRj  +  ipy,  where  pr^,  >  0,  and  at  v  =  ii»io,  at  which  pjo  may  be 
infinite.  Ordinarily  the  KK  relation  is  used  to  find  4»  from  the  measured  values  of  the  relative 
reflectance  R,  and  for  these  situations,  one  substitutes  Eq.  (46),  |  roi2 !  =  (R  Ro2)  w»  into  Eq. 
(72). 
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There  are  several  items  to  keep  in  mind  when  applying  Eq.  (72)  to  experimental  data.  If 
there  are  no  finite  zeros  of  ?oi2  on  the  positive  imaginary  axis,  the  substitution  =  00  gives 
-  it  for  -2  tan  - 1  yjo/v.  If  there  is  a  multiple  number  of  zeros  v0(,  where  vot- 1  is  less  than  t>0( 
along  the  positive  imaginary  axis,  the  summation  2  L(-  l)f  tan-1  vot  iv  should  replace 
-2  tan-1  vio/v. 

One  should  also  add  2v  to  the  right-hand  side  of  Eq.  (72)  for  every  ?Rj  that  is  less  than  v, 
if  the  principal  value  of  the  arctangent  function  is  used. 

3.1.8  Subtractive  Kramers-Kronig  (SKK)  Relation  for  4>(p) 

The  subtractive  Kramers-Kronig  (SKK)  relation  for  <f>  depends  upon  the  phase  0(vo) 
found  at  v  =  i-o  directly  from  an  experiment  or  computed  from  the  equation 

<M>o)  =  Im[fn  r0i2(v0)]  (73) 

after  one  enters  the  values  of  the  optical  constants  at  v0>  determined  in  a  separate 
experiment,  into  Eqs.  (4),  (8),  and  (35)(Ref.  19).  Usually  an  experimenter  chooses  vq  to  be  in 
a  transparent,  or  nearly  transparent,  region  somewhere  in  the  midst  of  the  reflectance 
spectrum.  As  noted  before,  Ahrenkiel  (Ref.  19)  has  shown  the  SKK  relation  to  be  less 
sensitive  than  the  corresponding  KK  equation  to  errors  introduced  by  the  extrapolations  of 
the  data  which  must  be  made  outside  the  data  range. 

3. 1.8.1  Derivation  of  the  SKK  Relation 

At  the  wavenumber  v$  at  which  <j>  is  known,  Eq.  (72)  is 

,/  v  a  ~2po  f  fn|roi2(vR)|  dvR 

<Kvo)  +  2m  t  =  - P  )  - - - - - 

IT  0  vft  -  V(5 

+  2  t  ( tan-i  -  tan-1  — 3i— ' ^ 

J"1  y  ^Rj  -  J'o  "Rj  +  "0  J 


-  2  tan-1 


no 


(74) 


Equations  (72)  and  (74)  can  be  combined  to  obtain  the  SKK  relation  for  <*>(v)  which  is 
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«->  +  2m.  =  ifcl  P  f  +  Z  [*(„„)  +  2mi] 

T  0  (4-  v2){v\  -  pj)  p0  L  J 


j  =  i 


[ 


+  2  .£.  I  tan-1— 2! - tan-1— 2L 


"Rj  -  V 


PRj  +  V 


-L  (tan 
»o  V 


-i _ v±  "U 


^Ij  -  *>0 


-  tan' 


"Rj  +  "o 


)] 


tan- 


V 

- tan-1 

PO 


(75) 


As  was  done  for  the  KK  relation,  Eq.  (72)  (R  Ro2)/l  can  be  substituted  for  'rm\.  All  the 
considerations  about  possible  zeros  on  the  positive  imaginary  axis  which  apply  to  the 
KK  relation,  also  apply  to  Eq.  (75).  They  appear  in  3. 1.7.5. 


3.2  LOCATION  OF  THE  ZEROS  OF  r012 


The  dispersion  relations  for  <t>,  the  phase  of  ron,  are  useless  unless  the  zero  points  of  r0|2 
are  known  in  the  upper  half  of  the  v-plane.  There  are  some  restrictions  on  the  form  that  the 
zeros,  Pj,  of  r0i2  can  have.  For  instance,  in  theory,  zeros  not  on  the  imaginary  axis  always 
appear  in  pairs  that  are  symmetric  in  position  relative  to  the  imaginary  axis  (Section  3.1 .7.3). 

In  this  section,  methods  are  examined  for  estimating  the  location  of  the  zeros  of  r0i2(p). 
The  general  scheme  is  to  estimate  the  functional  form  that  optical  quantities  have  on  the  real 
(pr)  axis  and  to  use  the  analytic  extension  of  these  forms  into  the  complex  v-plane  to 
compute  the  values  of  the  optical  quantities  away  from  the  real  axis. 

3.2.1  Finite  Number  of  Zeros 


Fortunately,  the  number  of  zeros  of  foi2  is  finite  because  of  the  forbidden  wavenumber 
gap  in  solid  dielectrics  and  semiconductors  (Ref.  48).  Near  the  gap  the  absorption  index  of 
the  film  greatly  exceeds  the  refractive  index,  and  the  zeros  of  r0i2  cross  into  the  lower  half  of 
the  complex  plane  in  the  manner  explained  by  Lupaskho  et  al.  (Ref.  7). 
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* 

3.2.2  Approximate  Values  of  i>i  and  n2 

Good  approximation  formulas  must  be  available  for  the  media  indices  ho,  r»i,  and  h2  if 
one  hopes  to  find  the  zeros  of  foi2-  The  assumption  is  made  that  medium  0  is  a  vacuum; 
hence,  ri0(»')  is  unity  for  all  complex  wavenumbers  v.  Thus,  approximate  functional  forms 
only  for  the  film  and  metal  indices,  hi  and  h2,  are  sought. 

In  the  preliminary  work,  the  analytic  extension  of  Eq.  (40), 

n2(i)  =  [1  -  r02  (*0]/[l  +  f»n(?>l  (76) 

was  used  to  compute  the  metal  index  from  the  extension  of  an  empirical  relation  for  r02,  such 
as  Eq.  (39), 

r<n(S)  =  -[0.939  -  (4.25  x  10" 5  cm);]*  (77) 

Note  that  Eq.  (77)  does  not  satisfy  the  relation  r02(  -  ?*)  =  roK*')  in  Section  3. 1 .7. 1 ,  and  the 
location  of  the  zeros  obtained  from  Eq.  (77)  are  not  paired  as  theory  predicts  in  Section 
3.1. 7.2. 

In  the  remainder  of  this  section  approximate  expressions  will  be  derived  for  the  optical 
constants  of  the  film,  n(j')  and  k(v).  First  consider  spectral  domains  where  the  film  causes 
noticeable  absorptions  of  radiation.  In  the  actual  wavenumber  domains  where  the  film  is 
sufficiently  absorbing,  the  film  acts  very  much  like  a  metal,  [  foi  |  is  much  greater  than 
|r12exp(2i71d1)  | ,  and  Eqs.  (4),  (34),  and  (35)  give 

Roi2  -  Inn  I2  =  [(1  -  n)2  +  k2]/[(l  +  n)2  +  k2]  (78) 

Solving  Eq.  (78)  for  k, 

k  «  {[Roi2(l  +  n)2  -  (1  -  n)2]/[l  -  R012]}*  (79) 

from  which  an  approximate  k-spectrum  can  be  constructed  if  crude  estimates  of  n  and 
experimental  values  of  R012,  are  obtained,  perhaps,  from  Eq.  (46)  using  measurements  of  the 
absolute  reflectance  R. 

The  estimates  of  the  k-  and  n-spectra  can  be  refined  if  Eqs.  (21)  and  (79)  are  incorporated 
in  an  iterative  process.  In  the  .first  iteration  n  may  be  set  to  unity  which  makes  Eq.  (79) 
become 
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k  *  2R^j2/(1  —  R012) (80) 

These  k-values  can  be  inserted  into  Eq.  (21)  to  find  an  n  -spectrum.  Subsequent  iterative 
steps  can  use  Eqs.  (21)  and  (79)  to  find  updated  k-  and  n-values. 

There  are  two  major  reasons  why  a  film  in  certain  spectral  regions  may  not  be 
sufficiently  absorbing  enough  to  be  represented  by  Eqs.  (78),  (79),  or  (80).  First,  the 
absorption  coefficient  k  of  the  film  medium  may  be  too  small  for  those  wavenumbers.  In 
those  spectral  ranges,  setting  k  to  zero  and  n  to  the  no-value  of  Eq.  (21)  should  give  good 
results.  The  second  reason  for  small  film  absorptions  is  that  the  film  thickness  d|  is  so  small 
that  |r12exp(2i71di)  is  comparable  to,  or  exceeds,  f  r0i  |  even  though  the  k-values  may  be 
large  for  those  wavenumbers.  The  construction  of  approximate  k-  and  n-spectra  is 
somewhat  more  difficult  in  those  spectral  domains  because  a  researcher  must  make  more 
assumptions  about  k  and  dj.  A  test  of  Eqs.  (79)  and  (80)  was  conducted  with  theoretical 
spectra,  and  results  are  discussed  in  Section  3.2.3. 

As  an  approximation  to  the  analytic  extension  of  the  index  of  the  film,  one  can  use 

(s’)  =  [ej(p)]  ^  (81) 

where  the  dielectric  constant  of  the  film  is  (Ref.  31) 

e,(v)  =  €0  +  Li  ap(*2  -  -  ibp?)-'  (g2) 

It  is  noted  that  e0  and  the  positive  constants  ap  and  bp  can  be  found  from  the  approximate 
k-and  n-spectra  of  Eqs.  (79)  and  (21),  by  computing  ^(v)  from 

ej  (?)  =  [n(v)  +  ik(v)]2  (83) 

then  fitting  ej(i^)  to  the  form  that  Eq.  (82)  has  along  the  real  axis,  i.e.,  when  v  =  v. 

3.2.3  Test  of  an  Iterative  Method  to  Find  the  Zeros 

Equation  (53)  has  been  incorporated  in  an  iterative  method  of  locating  the  zeros  of  roi2- 
For  each  value  of  the  index  j,  the  iterations  are  started  by  inserting  into  Eq.  (53)  the  roi,  r^, 
and  ni  values  computed  at  an  initial  .guess  vj(0)  of  the  zero,  thus,  getting  a  new  value. 
New  values  of  roi,  h2»  and  ni  at  pjd)  are  computed  to  reinsert  into  Eq.  (53)  for  the 
second  iteration.  The  iterations  are  repeated  until  the  moduli  of  successive  iterative  values 
differ  by  less  than  10  _  12  or  a  maximum  number  of  iterations  is  completed. 
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The  method  for  convergence  was  tested  using  theoretical  spectra  for  the  film  and  metal 
indices  n}  and  n2. 

A  theoretical  film  index  ni  was  calculated  from  Eqs.  (81)  and  (82).  The  choices  for  the 
parameters  of  Eq.  (82)  were  eo  =  1 .5  +  iO,  N  =  2  (the  number  of  absorption  bands),  ai 
=  a 2  =  3.0  x  10s  cm-2,  =  2000  cm-1,  v2  =  3000  cm-1,  and  bj  and  b2  =  10  cm-1.  It  is 

noted  that  n2(v)  was  obtained  from  Eq.  (76)  by  using  relations  similar  to  Eq.  (77)  for  ro2(»'). 

In  the  test  of  the  iterative  algorithm,  the  initial  guess  v(0)  was  set  to  be  the  origin,  film 
thickness  (dj)  values  varied  from  a  0.1  to  10  pm,  and  the  integer  j  progressed  over  the  values 
0,  ±1,  ±2,  until  the  converged  value  of  |  |  exceeded  10000  cm  -  *.  Calculations  were 

stopped  there  because  zeros  with  real  parts  that  are  far  from  the  data  wavenumber  range  give 
a  nearly  constant  contribution  to  the  phase  in  Eqs.  (72)  and  (75),  and  their  exact  locations 
would  not  be  important  in  a  typical  experiment. 

In  nearly  every  instance,  a  converged  value  of  'vx  was  found  which,  when  substituted  into 
the  analytic  extension  of  Eq.  (35),  yielded  a  value  of  roi2(»f)  whose  real  and  imaginary  parts 
differed  from  zero  by  less  than  10-6.  In  one  calculation,  j  was  0,  and  the  convergence  was 
not  as  good,  but  the  computed  zero  was  nearly  identical  to  that  for  j  =  - 1 .  In  all  of  the 
other  exceptional  cases,  the  j- values  were  negative  or  zero,  and  the  converged  zero  values 
were  in  the  second  quadrant  or  the  lower  half  of  the  complex  plane.  Zeros  in  the  first 
quandrant  did  not  seem  to  be  missing. 

Theoretically,  the  zeros  of  roi2  should  appear  in  pairs  reflected  about  the  imaginary  axis 
(Section  3. 1.7.2),  but,  as  was  pointed  out  in  Section  3.2.2,  the  use  of  an  empiricial  relation 
for  ro2»  such  as  Eq.  (77),  will  destroy  this  pairing.  In  fact,  when  Eq.  (77)  is  used  one  obtains 
an  approximate  pairing  of  the  zeros.  However,  when  r02(f')  =  -0.9,  so  that  r02(  -v*) 
=  roi2(*'*)»  the  converged  values  of  the  zeros  were  paired  according  to  =  -  v^\. 

The  same  hypothetical  spectra  were  also  used  to  compare  the  approximations  for  the  film 
absorption  coefficient  k  [Eqs.  (79)  and  (80)]  to  the  theoretical  values.  For  hypothetical  film 
thicknesses  1  jim  or  greater,  there  was  agreement  to  at  least  four  significant  digits  of  the 
theoretical  k-values  with  the  same  values  estimated  by  Eq.  (79)  in  the  strong  absorption 
regions.  Under  the  same  conditions,  the  values  found  by  Eq.  (80)  were  within  30  percent  of 
the  theoretical  values.  As  expected,  the  agreement  worsened  in  weaker  absorption  regions  or 
for  smaller  film  thicknesses  because  the  conditions  under  which  Eqs.  (79)  and  (80)  were  valid 
were  not  met  (Section  3.2.2). 
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3.3  COMPUTER  PROGRAM  SKKREFL 

The  authors  are  now  engaged  in  writing  and  testing  computer  program  SKKREFL  for 
the  SKK  analysis  of  experimental  reflectance  values  of  a  single-film  thickness  when  the  beam 
is  normally  incident  in  a  vacuum  upon  a  thin  film  deposited  on  a  metal  substrate.  The  three 
major  segments  of  the  program,  MAIN,  KAMKON,  and  NLNSYS,  are  essentially  the  same 
as  the  corresponding  segments  of  the  program  SKKTRANS  which  has  successfully  analyzed 
transmittance  data  for  a  single-film  thickness  (Section  2.2). 

KAMKON  does  the  numerical  integration  necessary  to  evaluate  the  SKK  dispersion 
relation  [Eq.  (75)]  for  the  phase  $  of  r0i2  at  a  given  wavenumber. 

NLNSYS  uses  the  experimental  value  of  r0l2  that  KAMKON  finds  from  Eqs.  (46)  and 
(47).  It  is  identical  to  the  segment  of  the  sarfle  name  in  SKKTRANS,  as  is  the  coding  needed 
to  choose  the  most  likely,  or  “best”  pair  of  n-  and  k-values  should  the  algorithms  yield  more 
than  one  pair  (Section  2.2). 

It  may  be  necessary  to  revise  the  film  optical  constants,  n  and  k,  obtained  from  NLNSYS 
because  the  algorithm  in  MAIN  which  locates  the  zero  points  of  r<>i2  uses  approximations  to 
these  constants.  The  program  can  use  the  n-  and  k-values  found  in  NLNSYS  in  a  search  for 
revised  zero  locations  and  compare  them  with  those  found  originally  in  MAIN.  If  a  revision 
is  necessary,  SKKREFL  can  redo,  one  or  more  times,  the  steps  in  KAMKON  and  NLNSYS 
with  updated  data  from  the  previous  pass. 

At  present,  program  SKKREFL  is  not  completely  assembled  or  fully  tested.  However, 
many  of  its  algorithms  have  worked  successfully  in  other  programs  previously  written,  and 
the  method  of  finding  zeros  of  roiz  has  passed  all  of  the  tests  to  date. 

4.0  THE  MULTIPLE-FILM  THICKNESS  METHOD;  PROGRAM  RENLIN 

Computer  program  RENLIN  was  written  to  analyze  the  normal  reflectance  of  a  beam 
incident  upon  a  film  deposited  on  a  metal  substrate  over  several  film  thicknesses.  RENLIN 
yields  the  complex  index  of  the  film  material,  nj,  which,  it  is  assumed,  does  not  depend  on 
film  thickness.  The  method  of  analysis  is  based  on  the  model  contained  in  Eqs.  (4),  (8),  and 
(34)  through  (44).  For  the  model,  it  is  noted  that  at  a  certain  wavenumber,  the  film  indices  n 
and  k,  and  all  other  parameters,  are  constant  by  assumption  during  an  experiment. 
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RENLIN  fits  the  experimental  reflectance  values  of  two  or  more  film  thicknesses  at  a 
given  wavenumber  to  the  model  of  Eqs.  (4),  (8),  and  (34)  through  (44).  The  fitting  procedure 
is  the  nonlinear  least-squares  algorithm  of  Marquardt  (Ref.  8),  and  the  fitted  parameters  are 
n  and  k  of  the  film  material.  Except  for  the  equations  representing  the  model  and  some 
modifications  of  the  data  input  procedures,  RENLIN  is  identical  in  coding  to  TRNLIN,  the 
program  that  fits  experimental  transmittance  values  of  thin  films  used  in  Ref.  1.  It  is 
recommended  that  the  present  version  of  the  nonlinear  least-squares  algorithm  in  program 
segment  NLIN,  which  is  a  part  of  both  of  the  RENLIN  and  TRNLIN  programs,  be  revised. 
The  coding  of  this  algorithm  is  in  single-precision  arithmetic,  whereas  all  the  other 
algorithms  use  double  precision.  It  also  happens  to  be  very  confusing  to  follow.  The  main 
reason  for  the  proposed  revision  of  NLIN,  however,  is  to  improve  jts  efficiency.  Least- 
squares  algorithms  tend  to  be  costly  procedures  involving  minutes  of  CPU  time  and 
kilobytes  of  memory.  Huber  (Ref.  49)  has  recently  discovered  a  new  “extrapolated  least- 
squares  optimization”  method  that  typically  reduces  the  number  of  iterative  steps  to  about 
half  that  of  the  conventional  algorithms.  The  method  is  actually  a  slight  modification  of  the 
present  algorithm  in  NLIN  and  would  improve  each  of  the  present  iterative  steps. 

5.0  CONCLUSIONS 

• 

The  lifetime  of  space  satellites  is  often  directly  related  to  contamination  of  critical  optical 
surfaces  by  condensed  gas  films.  A  basic  property  required  for  determining  optical 
contamination  effects  is  the  complex  refractive  index  of  the  gas  film.  Results  of  efforts  to 
develop  mathematical  models  and  computer  programs  for  the  determination  of  the  complex 
mathematical  models  and  computer  programs  for  the  determination  of  the  complex 
refractive  index  from  optical  transmittance  and  reflectance  measurements  are  presented  in 
this  report. 

In  previous  work  (Ref.  1),  the  program  TRNLIN  has  been  used  successfully  to  determine 
refractive  index  values  using  experimental  transmittance  data  for  multiple-film  thicknesses 
of  condensed  gases.  The  program  RENLIN  has  been  written  to  analyze  the  normal 
reflectance  of  a  beam  incident  on  a  film  deposited  on  a  metal  substrate  over  several  film 
thicknesses.  RENLIN  yields  the  complex  index  of  the  film  material  which  is  assumed  not  to 
depend  on  film  thickness.  Except  for  equations  representing  the  reflectance  model  and  some 
modifications  of  the  data  input  procedures,  RENLIN  is  identical  in  coding  to  TRNLIN. 

Detailed  descriptions  of  two  new  programs  to  determine  the  complex  refractive  index 
from  experimental  measurements  of  a  single-film  thickness  are  presented.  The  advantage  of 
this  model  is  that  the  assumption  used  in  previous  models  that  the  optical  properties  of  a 
dielectric  film  are  independent  of  film  thickness  is  not  required.  The  program  SKKTRANS 
utilizes  a  new  subtractive  Kramers-Kronig  algorithm  which  gives  highly  reliable  index  values 
from  measurements  of  the  absolute  transmittance  of  a  single  uniform  thin  film  deposited  on 
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a  uniform  thick  substrate.  This  program  reduces  the  number  of  necessary  experimental 
measurements  and  is  capable  of  noting  any  differences  in  the  optical  constants  of  the  film 
material  due  to  changes  in  film  thickness.  A  theoretical  CO^-cryofilm  transmittance 
spectrum  was  used  to  test  the  algorithms  in  SKKTRANS.  Excellent  agreement  was  found 
between  calculated  and  original  index  values.  From  the  results  of  the  SKKTRANS  program 
and  the  transmittance  data  taken  previously  from  an  N2/CO2  mixture,  it  was  observed  that 
the  absorption  index  determined  through  the  SKKTRANS  program  increased  with 
thickness.  Further  investigation  is  required  to  determine  whether  this  trend  is  a  real  physical 
effect  or  a  computational  artifact.  This  same  effect  was  observed  for  other  gas  mixtures.  The 
values  computed  for  individual  thickness  using  SKKTRANS  were  in  close  agreement  with 
those  obtained  with  TRNLIN  and  multiple-film  thicknesses.  Also  the  program  SKKREFL 
was  derived  for  determining  n  and  k  using  subtractive  Kramers-Kronig  technique  from 
experimental  reflectance  data  for  a  single-film  thickness.  This  program  is  not  completely 
assembled  and  has  not  been  fully  tested. 
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